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ABSTRACT 

The  purpose  of  this  report  is  to  present  the  additional  results  of 
applying  the  Point  Scattering  Model  for  radar  image  simulation  to  three 
(3)  new  cases.  This  work  is  supplemental  to  that  previously  reported  ’ . 
The  work  was  sponsored  by  the  1).  S.  Army  Enqineer  Topoqraphic  Laboratories 
(ETL). 

The  results  reported  were  obtained  for  four  (4)  simulations  corres- 
ponding to  four  specific  altitudes  in  the  terminal  trajectory  of  a guided 
missile  (three  of  the  simulations  are  new),  each  successively  lower.  The 
sequence  of  four  simulations  was  tested  against  actual  radar  data  of  the 
same  site  via  the  Correlatron*.  The  correlation  tests  showed  performance 
of  simulations  produced  by  the  PSM  to  be  acceptable  in  three  cases,  and  to 
be  unacceptable  in  one  case.  The  unacceptable  performance  was  explained 
by  geometric  errors  inadvertently  introduced  into  the  data  base  and  was  not 
produced  by  the  simulation  method. 
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This  document  was  prepared  by  the  Kansas  Simulation  Group,  Remote 
Sensing  Laboratory  (RSL),  The  University  of  Kansas,  Lawrence,  Kansas,  to 
report  the  work  performed  and  results  obtained  for  a radar  simulation  study 
performed  under  Contract  DAAG53-76-C-0154,  modification  P00003,  dated 
22  June  1977,  with  the  Engineer  Topographic  Laboratories  (ETL),  Fort  Belvoir, 
Virginia. 

The  original  effort  performed  was  a radar  simulation  study  to  validate 
the  Point  Scattering  Method,  a radar  image  simulation  model  which  had  been 
developed  previously,  to  investigate  terrain  feature  extraction  techniques 
for  constructing  category  data  bases  (digital  ground  model)  for  radar 
simulation,  and  to  use  this  point  scattering  radar  simulation  model  to 
generate  radar  reference  scenes  for  terminal  guidance  applications.  The 
work  performed  and  results  obtained  in  the  original  study  have  been  reported 
previously  ^ . 

The  study  conducted  under  modification  P00003  included  preparation  of 
three  (3)  radar  reference  scenes  at  three  (3)  different  scales.  The  work 
performed  to  construct  these  three  additional  radar  reference  scenes  and 
results  obtained  from  using  them  in  a test  configuration  are  reported  in 
this  document. 

This  document  is  designed  to  be  an  addendum  to  the  reports  prepared 
for  the  original  effort1’^  and  reports  only  the  additional  work  and  results 
obtained  under  modification  P00003.  It  reports  only  the  additional  results 
and  relies  heavily  upon  the  knowledge  of  the  earlier  documents  for  philo- 
sophical discussions,  derivations,  and  technical  details. 
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PURPOSE 


The  purpose  of  this  document  is  to  report  the  results  of  applying 
the  Point  Scattering  Method  (PSM),  a radar  image  simulation  model,  to 
construction  of  radar  reference  scenes  for  a missile  terminal  guidance 
application  utilizing  area  correlation  as  the  guidance  technique.  The 
work  was  performed  at  the  Remote  Sensing  Laboratory  (RSL),  The  University 
of  Kansas,  Lawrence,  Kansas  under  contract  with  the  Engineer  Topographic 
Laboratories  (ETL),  United  States  Army,  Fort  Belvoir,  Virginia. 

The  purpose  of  the  work  performed  was  to  test  the  feasibility  of 
using  the  PSM,  previously  developed  and  reported  ’ , to  produce  high 
quality  radar  reference  scenes  for  terminal  guidance  applications.  This 
purpose  was  accomplished  by  constructing  the  three  (3)  additional  radar 

1 2 

reference  scenes  needed  to  go  with  the  one  (1)  developed  in  previous  work 
to  make  a complete  set  for  one  (1)  target  site. 
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The  scope  of  the  work  performed  in  this  feasibility  study  regarding 
the  use  of  the  PSM  to  make  high  quality  radar  reference  scenes  for  terminal 
guidance  applications  was  limited  to  constructing  three  (3)  radar  reference 

» 

scenes  for  one  (1)  target  site.  The  radar  reference  scene  previously 

1 2 

constructed  and  reported  ’ was  the  largest  scale  (and,  thus,  greatest 
detail)  scene  required  in  a set  of  four  (4).  The  three  (3)  radar  refer- 
ence scenes  produced  in  the  work  reported  here  were  the  successively  smal- 
ler scale  (and,  thus,  larger  area  coverage)  scenes  in  the  set. 

The  lone  target  site  selected  for  the  feasibility  study  of  the  PSM 
was  the  topographic  region  in  the  states  of  Alabama,  Tennessee,  and 
Mississippi  centered  on  the  northwest  corner  of  the  power  house  at  the  Pick- 
wick Landing  Dam,  Tennessee  (coordinates;  34°  04'  15"  N by  88°  15'  05"  W). 


EXECUTIVE  SUMMARY 

The  Point  Scattering  Model  (PSM)  for  radar  simulation,  developed  at 
the  RSL  (Remote  Sensing  Laboratory,  University  of  Kansas,  Lawrence, 

Kansas)  has  been  applied  to  the  problem  of  synthesizing  radar  reference 
images  for  use  in  the  terminal  guidance  system  of  a missile.  The  guidance 
system  for  which  the  PSM  was  tested  utilizes  a Correlatron*  for  deriving 
guidance  information  from  the  two-dimensional  cross  correlation  between 
reference  images  and  actual,  "live"  radar  data. 

Four  (4)  reference  images  were  produced  via  the  PSM,  each  successively 

representing  lower  altitudes  and  later  stages  of  terminal  descent  over  a 

target.  The  two  higher  altitude  reference  images  were  produced  via  the 

PSM  from  one  radar  data  base,  and  the  two  lower  altitude  reference  images 

from  a second  data  base.  The  first  of  these  data  bases  was  constructed 

in  work  performed  and  reported  in  this  document.  The  second  data  base 

1 2 

was  constructed  in  earlier  work  ’ for  band  1 simulations,  and  was  upgraded 
here  for  both  bands  1 and  2 simulations.  The  same  philosophy  and  construc- 
tion techniques  were  utilized  for  both  data  bases.  The  bands  3 and  4 data 
base  covers  a much  larger  geographic  area  and  was  built  having  a smaller 
inherent  scale  than  the  bands  1 and  2 data  base,  consistent  with  the 
higher  altitude  simulations  it  is  designed  to  support.  Aside  from  these 
obvious  differences  in  area  or  coverage  and  scale,  a major  difference  exists 
between  the  two  data  bases.  The  bands  3 and  4 data  base  was  built  using 
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a 1:100,000  scale  ortho  photo  and  elevation  data  developed  via  UNAMACE* 
as  the  input  source  intelligence  information.  The  bands  1 and  2 data 
base  was  constructed  using  standard  7-1/2  quadrangle  USGS  (United  States 
Geological  Survey)  1:24,000  scale  maps  and  elevation  data  derived  from 
the  maps  by  DMA  (Defense  Mapping  Agency)  as  the  input  source  intelligence 
data.  The  ortho  photo  was  constructed  as  a tangent  plane  projection 
centered  on  the  target  site.  The  USGS  maps  are  modified  polyconic  projec- 
tions. Thus,  differences  in  geometric  fidelity  underlie  the  two  data  bases. 

The  sequence  of  four  (4)  reference  scenes  was  tested  on  a Correlatron 
at  ETL  (U.S.  Army  Engineer  Topographic  Laboratories)  versus  actual  radar 
data  previously  collected  and  recorded  in  a flight  test  program  conducted 
over  the  Pickwick  test  site.  Results  of  this  test  were  very  favorable 
for  the  PSM,  but  results  were  not  as  decisive  as  expected  because  of 
geometric  fidelity  problems  inadvertently  created  by  mechanical  means 
when  constructing  the  bands  3 and  4 data  base.  A qualitative  assessment 
of  the  results  produced  is: 

1)  Band  1 (Lowest  Altitude)  - Very  good  performance 
within  accuracy  requirements; 

2)  Band  2 (Higher  Altitude)  - Excellent  performance 
well  within  desired  accuracies; 

3)  Band  3 (Next  Higher  Altitude)  - Failure  to  perform 
within  accuracy  requirements; 

4)  Band  4 (Highest  Altitude)  - Acceptable  performance 
within  accuracy  requirements. 

These  results  represent  a single  preparation  of  the  radar  reference 
images  and  are  not  the  best  possible  results  produced  after  much  changing 
and  custom-tailoring  of  various  parameters  as  this  has  not  been  done. 

The  bands  1 and  2 simulations  were  developed  from  the  data  base  having 
the  highest  degree  of  detail.  Performance  of  these  simulations  was 
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spectacular  in  view  of  the  fact  that  they  worked  so  well  on  the  first 
trial.  Band  3 and  4 simulations  were  developed  from  the  data  base  having 
geometrical  problems  as  discussed  in  Section  6.  Performance  of  the 
simulations  produced  from  this  data  base  was  better  than  expected  in  view 
of  the  nature  of  the  problems  built  into  the  data  base. 

The  problems  in  the  bands  3 and  4 data  base  are  unfortunate  because 
they  cloud  what  otherwise  should  have  been  a clear  decision  in  favor  of 
the  PSM.  Surely  if  this  data  base  had  not  contained  geometrical  errors, 
the  PSM  simulations  for  bands  3 and  4 would  have  performed  as  well  as 
the  bands  1 and  2 simulations  did.  Of  course,  this  cannot  be  proved  by 
the  data  that  exist,  but  it  is  a reasonable  conjecture,  especially  in 
view  of  the  fact  that  even  with  the  errors,  the  simulation  for  band  4 
was  acceptable. 

Simulations  produced  via  the  PSM  should  work  well  in  this  guidance 
application.  The  PSM  rigorously  models  the  radar  guidance  system,  and 
accuracy  of  performance  in  this  application  should  be  limited  by  the 
accuracy  of  the  data  base. 

In  conclusion,  the  results  produced  in  work  performed  in  this  study 
verify  the  validity  of  using  the  PSM  for  producing  radar  reference  scenes 
to  be  used  in  conjunction  wi th  th?  Correlatron.  The  verification  is  not 
as  strong  as  it  might  have  been  if  errors  of  construction  technique  (i.e., 
implementation  or  mechanical  errors,  not  errors  in  construction  philosophy) 
had  not  been  introduced  into  the  bands  3 and  4 data  base.  The  work 
performed  show  the  PSM  to  be  a viable  tool  for  solving  "real-world"  prob- 
lems, not  just  an  academic  exercise,  which  can  and  should  be  applied  to 
many  present  and  future  problems. 
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1.0  INTRODUCTION 

The  work  reported  in  this  document  was  performed  to  test  the 
feasibility  of  using  the  Point  Scattering  Method  (PSM),  a previously 
developed,  validated,  and  reported  radar  image  simulation  model  J’2,3’^ 
to  produce  high  quality  radar  reference  scenes  for  terminal  guidance 
applications.  The  study  performed  here  was  supplemental  to  the  work 
previously  performed  and  reported  in  references  3 and  4.  In  the  pre- 
vious work  preliminary  indications  of  using  PSM,  for  terminal  guidance 
applications  were  investigated  by  producing  one  (1)  radar  reference  scene 
for  one  (1)  target  site.  That  test  was  very  successful  but  was  deemed  not 
conclusive  because  the  terminal  guidance  application  for  which  PSM  was 
being  tested  required  a set  of  four  radar  reference  scenes  for  each 
target  site  and  only  one  (1)  of  the  set  (the  lowest  altitude  and  largest 
scale)  had  been  produced  by  the  PSM.  The  work  reported  here  was  performed 
to  produce  the  additional  three  (3)  radar  reference  scenes  required  to 
complete  the  set  for  the  one  target  site.  The  target  site  for  which 
this  sequence  of  reference  scenes  was  formed  was  the  Pickwick  Landing 
Dam  test  site.  This  site  is  located  in  the  states  of  Alabama,  Tennessee, 
and  Mississippi  with  the  site  centered  on  the  northwest  corner  of  the  power 
house  at  the  dam  (coordinates  34°  04'  15"N  by  88°  15'  05"  W). 


^Holtzman,  J.  C.,  V.  H.  Kaupp,  J.  L.  Abbott,  V.  S.  Frost,  E.  E.  Komp,  and 
E.  C.  Davision,  "Radar  Image  Simulation:  Validation  of  the  Point  Scattering 
Method,"  Vol.  I,  ETL-0117,  The  University  of  Kansas  Center  for  Research, 
Inc.,  September,  1977,  AD-A053  253. 

2Holtzman,  J.  C.,  V.  H.  Kaupp,  J.  L.  Abbott,  V.  S.  Frost,  E.  E.  Komp,  and 
E.  C.  Davison,  "Radar  Image  Simulation:  Validation  of  the  Point  Scattering 
Method,"  Vol.  II,  ETL-0018  , The  University  of  Kansas  Center  for  Research, 
Inc.,  September,  1977,  AD-A053  240. 

3Holtzman,  J.  C.,  V.  H.  Kaupp,  R.  L.  Martin,  E.  E.  Komp,  V.  S.  Frost, 

"Radar  Image  Simulation  Project:  Development  of  a General  Simulation  Model 
and  an  Interactive  Simulation  Model,  and  Sample  Results,"  ETL-0047,  RSL 
Technical  Report  234-13,  The  University  of  Kansas  Center  for  Research, 

Inc.,  February,  1976,  AD-A027  151. 
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Holtzman,  J.  C.,  V.  H.  Kaupp,  and  J.  L.  Abbott,  "Radar  Image  Simulation 
Project,"  ETL-0098,  RSL  Tech.  Report  234-15,  University  of  Kansas  Center 
for  Research,  Inc.,  September,  1976,  AD-A051  501. 


Philosophical  discussions  of  the  terminal  guidance  problem,  the  PSM 
radar  image  model,  and  work  plan  are  presented  in  this  section  (1.1,  1.2, 
and  1.3,  respectively)  and  the  technical  details  of  the  work  performed 
and  the  results  obtained  are  incorporated  in  succeeding  sections. 


1.1  Terminal  Guidance 


< 
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The  task  of  delivering  an  explosive  projectile  accurately  enough 
to  destroy  an  intended  target  has  three  (3)  basic  components:  (1)  Distin- 
guish the  target  from  its  surroundings;  (2)  Determine  the  exact  position 
of  the  target  relative  to  the  launching  point  of  the  missile;  (3)  Guide 
the  projectile  onto  the  target.  "Terminal  guidance"  is  the  descriptive 
phrase  used  to  signify  the  final  segment  of  an  explosive  projectile's 
flight  in  which  it  is  guided  onto  the  target  (item  3,  above).  In  this 
work  we  are  concerned  only  with  the  question  of  "terminal  guidance". 

How  the  first  two  (2)  phases  of  delivering  a projectile  onto  a target 
(items  1 and  2,  above)  are  satisfied  doesn't  enter  thi *-  concern.  A 
variety  of  ways  exist  to  deliver  a projectile  into  the  neighborhood  of 
a target.  These  span  the  gamut  from  ballistic  to  cruise  to  bomber  delivery 
techniques. 

Different  terminal  guidance  philosophies  exist.  The  philosophy 
for  which  the  feasibility  of  using  the  PSM  was  investigated  was  area 
correlation.  In  this  particular  scheme  reference  radar  images  of  the 
target  site  (computer- generated  simulations  formed  via  the  PSM  in  this 
work)  are  stored  on-board  the  projectile  in  its  guidance  package.  In 
the  terminal  phase  of  its  trajectory  a radar  interrogates  the  terrain 
below  the  projectile  and  the  guidance  device  compares  the  radar  data 
to  the  previously  stored  reference  scene, deriving  guidance  data  from 
the  comparison.  The  comparator  employed  in  the  work  performed  here 
was  the  Correlatron,  an  area-cross-correlation  measuring  device  manufac- 
tured by  Goodyear  Aerospace.  The  Correlatron  has  been  discussed  in  Avia- 

5 

tion  Week  and  in  the  earlier  report  of  applying  the  PSM  to  the  question 
of  terminal  guidance,6  thus  a discussion  of  the  Correlatron  will  not  be 
repeated  here. 


5 

Klass,  Philip,  J.,  "Guidance  Device  Set  for  Pershing  Tests,"  Aviation 
Week  and  Space  Technology,  May  12,  1975. 

6Holtzman,  J.  C.,  V.  H.  Kaupp,  J.  L.  Abbott,  V.  S.  Frost,  E.  E.  Komp,  and 
E.  C.  Davison,  "Radar  Image  Simulation:  Validation  of  the  Point  Scattering 
Method,"  Vol.  I,  ETL-0117,  The  University  of  Kansas  Center  for  Research, 
Inc.,  September,  1977,  AD-A053  253,  pp.  105-110. 
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The  guidance  package  incorporating  a Correlatron  requires  that  four 
(4)  reference  radar  images  be  stored  on-board  the  projectile  for  comparison 
with  “live"  video  radar  data  (acquired  via  a resident  radar)  of  the 
terrain  in  the  vicinity  of  the  intended  target  (this  presupposes  that 
the  first  two  phases  of  the  projectile  problem  have  been  satisfied 
accurately  and  the  projectile  is  in  the  neighborhood  of  the  target). 

These  four  reference  scenes  form  a sequence  starting  with  the  smallest 
scale  (largest  area  coverage)  and  progressing  to  the  largest  scale 
(smallest  area  coveraqe)  corresponding  to  the  highest  altitude  to  the 
lowest  altitude,  respectively. 

It  is  this  sequence  of  four  (4)  reference  scenes  for  which  the 
feasibility  of  using  the  PSM  was  tested  in  the  work  reported  here  and 
in  the  earlier  work^.  Radar  reference  images  have  been  formed  via 
the  PSM  in  this  work  for  the  three  (3)  smallest  scale  (each  progressive- 
ly larger)  scenes  corresponding  to  the  three  (3)  hiqhest  altitudes, 
and  a radar  reference  image  was  produced  in  the  earlier  work  via  the  PSM 
for  the  largest  scale  (lowest  altitude)  scene.  These  four  (4)  radar 
reference  images  toqether  comprise  a complete  sequence  of  scenes  for 
one  target  site. 


Holtzman,  J.  C.,  V.  H.  Kaupp,  J.  L.  Abbott,  V.  S.  Frost,  E.  E.  Komp,  and 
E.  C.  Davison,  "Radar  Image  Simulation:  Validation  of  the  Point  Scatter- 
ing Method,"  Vol.  I,  ETL-0117,  The  University  of  Kansas  Center  for 
Research,  Inc.,  September,  1977,  AD-A053  253. 
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Holtzman,  J.  C.,  V.  H.  Kaupp,  J.  L.  Abbott,  V.  S.  Frost,  E.  E.  Komp,  and 
E.  C.  Davison,  "Radar  Image  Simulation:  Validation  of  the  Point  Scatter- 
ing Method,"  Vol.  II,  ETL-0118,  The  University  of  Kansas  Center  for 
Research,  Inc.,  September,  1977,  AD-A053  240. 


1.2  The  Point  Scattering  Method:  A Model  for  Radar  Image  Simulation 


The  Point  Scattering  Method  (PSM)  for  radar  image  simulation  has  been 
3 1 2 

developed  and  validated  ’ in  previous  work.  A detailed  description  is 
not  repeated  here  but  the  essence  of  the  PSM  is  summarized  to  lay  the  frame- 
work for  understanding  what  is  involved  in  constructing  radar  reference 
scenes  for  the  Correlatron  guidance  package.  The  appropriate  references 
should  be  consulted  for  additional  technical  information. 


1.2.1 


Introduction  to  the  PSM 


The  PSM  is  a model  developed  for  the  end-to-end  simulation  of  radar 
images.  By  radar  image  simulation  is  meant  synthesis  of  the  image  which 
would  have  been  produced  by  a radar  had  it  been  flown  according  to  a 
prescribed  flight  path  over  the  actual  ground  site.  The  radar  senses  the 
terrain  in  the  microwave  portion  of  the  electromagnetic  spectrum  and  pro- 
duces an  image  on  ordinary  photographic  film.  The  radar  image  represents 
the  reflectivity  characteristics  of  objects  on  the  ground  at  radar  wave- 
lengths and  displays  the  terrain  in  fine  detail  and  with  spectacular  relief. 
So  does  the  simulated  radar  image. 


i 

| 
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Holtzman,  J.  C.,  V.  H.  Kaupp,  R.  L.  Martin,  E.  E.  Komp,  V.  S.  Frost, 

"Radar  Image  Simulation  Project:  Development  of  a General  Simulation 
Model  and  an  Interactive  Simulation  Model,  and  Sample  Results,"  ETL-0047, 
RSL  Technical  Report  234-13,  The  University  of  Kansas  Center  for  Research, 
Inc.,  February,  1976,  AD-A027  151. 

^Holtzman,  J.  C.,  V.  H.  Kaupp,  J.  L.  Abbott,  V.  S.  Frost,  E.  E.  Komp,  and 
E.  C.  Davison,  "Radar  Image  Simulation:  Validation  of  the  Point  Scattering 
Method,"  Vol.  I,  ETL-0117,  The  University  of  Kansas  Center  for  Research, 
Inc.,  September,  1977  , AD-A053-253. 
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Holtzman,  J.  C.,  V.  H.  Kaupp,  J.  L.  Abbott,  V.  S.  Frost,  E.  E.  Komp,  and 
E.  C.  Davison,  "Radar  Imaqe  Simulation:  Validation  of  the  Point  Scattering 
Method,"  Vol.  II,  ETL-01 18,  The  University  of  Kansas  Center  for  Research, 
Inc.,  September,  1977,  AD-A053  240. 
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The  Point  Scattering  Method  represents  the  radar,  the  terrain  and 
its  microwave  response,  and  the  image  and  its  record  on  photographic  film 
by  a closed-system  model  in  this  guidance  specialization  of  the  model. 

The  model  rests  on  firm  theoretical  foundations  and  is  mathematically 
rigorous.  It  incorporates  all  aspects  of  the  imaging  radar  problem.  The 
model  is,  however,  more  than  equations.  It  is  a philosophy  of  radar  image 
synthesis. 

Just  as  all  mathematical  models  are  abstractions  of  reality,  so  is 
the  PSM.  It  attempts  to  describe  in  closed  form  the  processes  of  the 
imaging  radar  system  consisting  of  the  radar,  ground,  and  image.  The 
PSM  is  completely  general  and  is  capable  of  producing  images  having 
whatever  accuracy  is  required  (at  least  for  distributed  targets)  and 
producing  them  in  a desired  output  format  for  any  particular  application 
if  the  cost  is  paid  in  time,  complexity,  and  resources  expended.  The 
true  value  of  the  PSM  comes  from  the  flexibility  to  specialize  it  for 
each  application  by  making  valid  approximations  and  simplifications.  The 
model  can  be  made  efficient  and  cost-effective  to  use. 

The  development  of  the  general  PSM  and  the  philosophy  of  radar  image 
simulation  embodied  by  this  model  are  summarized  according  to  the  follow- 
ing four  natural  divisions  of  the  problem:  (1)  Imaging  model;  (2)  Radar 
geometrical/propagation  phenomena;  (3)  Ground  model;  and  (4)  Reflectivity 
model.  In  all  that  follows,  it  should  be  kept  in  mind  that  regardless 
of  how  the  model  is  implemented  (e.g.,  on  a digital  computer),  there  are 
three  basic  input  data  requirements  of  the  model.  These  are  illustrated 
in  the  block  diagram  of  Figure  1.  First,  the  operating  parameters  of 
the  radar  being  simulated  must  be  incorporated.  Second,  a symbolic  repre- 
sentation of  the  relief  and  dielectric  properties  of  each  selected  ground 
site  must  be  specified.  This  symbolic  representation  of  the  terrain  (the 
ground  truth)  is  called  the  data  base.  Last,  the  reflectivity  properties 
of  all  the  different  dielectric  categories  (radar  categories)  included 
in  the  ground  truth  data  base  must  be  specified.  The  reflectivity  proper- 
ties of  the  different  categories  are  modelled  by  the  differential  scatter- 
ing cross-section  (o°);  either  theoretical  or  empirical  backscatter  data 
are  used. 


Upon  satisfaction  of  all  three  input  data  requirements,  the 
computational  algorithms  of  the  Point  Scattering  Model  are  employed  to 
synthesize  the  radar  image  that  would  have  been  collected  if  the 
modeled  radar  had  been  flown  over  the  selected  ground  site,  a fact 
attested  to  by  the  images  presented  later. 

Previous  efforts  in  the  area  of  radar  image  simulation  have 
incorporated  over-simplifications  in  either  the  model  of  the  imaging 
process  or  in  the  representation  of  the  microwave  reflectivity  of  the 
various  categories.  The  Point  Scattering  Method  does  not  rely  on 
arbitrary  over-simplifications;  it  is  completely  rigorous  and  general. 
Thus,  for  the  first  time  a complete  radar  image  simulation  model, 
representative  of  all  facets  of  the  process,  is  available. 


1.2.2  Imaging  Model 


The  imaging  model  is  the  final  computational  algorithm  of  the  Point 
Scattering  Method.  After  the  ground  truth  data  base  (terrain  feature 
model)  of  the  desired  site  has  been  specified,  the  reflectivity  (a°)  data 
for  the  various  categories  included  in  the  data  base  have  been  obtained, 
and  the  complex  geometry  relating  the  radar  platform  to  the  scene  has  been 
determined,  the  imaging  model  is  used  to  calculate  the  power  reradiated 
from  the  ground  back  to  the  radar  for  each  pixel  (picture  element)  in  the 
image.  It  calculates  the  intensity  of  the  ground-return  signal  exiting 
the  receiver,  and  in  the  present  guidance  specialization  of  the  PSM,  it 
calculates  the  conversion  of  this  intensity  into  density  of  silver  grains 
in  the  exposed  and  developed  image.  This  algorithm  determines  the  shade 
of  grey  of  each  pixel  in  the  image,  and  hence,  is  called  the  greytone 
equation.  The  greytone  equation  produces  the  final  results,  drawing  upon 
all  preceding  data  and  calculations.  It  relates  the  around  to  the  radar 
and  to  the  image. 

An  appropriate  starting  point  for  this  develonment  of  the  imaging 
model  would  be  a description  of  the  principles  of  imaging  radars.  However, 
as  these  principles  have  been  well  documented  in  the  literature,  they  will 
not  be  repeated  here^’®.  The  starting  point  for  the  radar  simulation 
imaging  model  (greytone  equation)  is  the  prediction  (estimation)  of 
the  power  reradiated  from  each  resolvable  ground  element  (resolution  cell) 
The  size  of  a resolution  cell  is  determined  in  the  range  direction  by  the 
pulse  length  and  in  the  azimuth  direction  by  the  directional  properties  of 
the  antenna,  either  real  or  synthetic.  For  the  purposes  of  radar  image 
simulation,  it's  assumed  that  the  ground  can  be  modeled  as  a collection  of 
homogeneous  regions,  each  at  least  the  size  of  a resolution  cell  (these  are 
called  distributed  targets).  Point,  or  cultural,  targets  can  be  handled 


7Skolnik,  M.  K. , (Editor),  Radar  Handbook.  New  York:  McGraw-Hill,  1970. 
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Moore,  R.  K.,  "Microwave  Remote  Sensors,"  in  R.  G.  Reeves,  Remote 
Sensing  Manual,  Falls  Church,  Virginia:  American  Society  of  Photo- 
grammetry,  1975,  Chapter  9. 
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statistically,  symbolically,  by  scattering  calculations,  or  alternately 

g 

by  the  area  spatial  filtering  method  . The  power  reradiated  from  all 
the  scattering  centers  located  within  a resolution  cell  combine  at  the 
receiving  antenna  to  produce  one  value  for  the  resolution  cell  which  is 
designated  the  return  power  (PR). 

In  reality  the  return  power  is  not  a deterministic  process  as  the 
foregoing  discussion  might  imply.  The  amplitude  of  the  return  power 
received  by  the  antenna  mounted  on  a moving  vehicle  fluctuates  widely 
because  of  variations  in  the  phase  of  the  reradiation  from  the  different 
scatterers  in  the  illuminated  area  (resolution  cell).  This  phenomenon 
accounts  for  the  speckled  nature  ("grainy"  appearance)  of  radar  images 
and  is  called  "fading".  The  statistics  of  a "fading"  signal  have  been 
well-documented  for  most  homogeneous  targets  and  it  has  been  shown  that 
the  signal  amplitude  can  be  described  by  a Rayleigh*  probability  distri- 


bution 


7,8,10 


^Frost,  V.  S. , J.  L.  Abbott,  V.  H.  Kaupp,  and  J.  C.  Holtzman,  "A  Mathe- 
matical Model  for  Terrain-Imaging  Radar  and  Its  Potential  Applications 
to  Radar  Imane  Simulation,"  University  of  Kansas  Center  for  Research, 
Inc.,  Remote  Sensing  Lab.,  Tech.  Report  319-6,  Lawrence,  Ks.,  Nov.  1976. 

*The  Rayleigh  probability  distribution  is  one  of  the  probability  models^’*’* 
which  can  be  used  to  describe  the  signal  amplitude  variation. 

Use  of  a different  probability  model  would  modify  the  resulting 
"fading"  per  cell  somewhat.  However,  for  medium  to  coarse 
resolution  radars,  final  averaging  yields  similar  results. 

7Skolnik,  M.  I.,  (Editor),  Radar  Handbook.  New  York:  McGraw-Hill,  1970. 
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Manual,  Falls  Church,  Virginia:  American  Society  of  Photogrammetry , 
1975,  Chapter  9. 
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Backscatter  from  Selected  Agricultural  Targets,"  IEEE  Trans.  Geosci. 
Electron.,  Vol . GE-13,  October  1976,  pp.  149-157. 
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cessors," J.  Opt.  Soc.  Amer. , Vol.  66,  No.  11,  November  1976. 
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If  square-law  detection  is  assumed  for  the  radar  being  modeled, 
then  the  postdetection  signal  is  a random  variable  having  a chi- 
square  probability  density  function  with  2NS  degrees  of  freedom1*1 
where  NS  specifies  the  number  of  "independent  samples"  being  averaged. 
The  minimum  width  of  a backscatter  lobe  in  the  azimuth  direction  (best 
azimuth  resolution)  is  specified  as  L/2  where  L is  the  real  antenna 
length11.  A real-aperture  SLAR  (Side-Looking  Airborne-Radar)  has  an 
azimuth  resolution  given  by  bR  where  B is  the  diffraction-limited 
beamwidth  given  by  x/L  (the  illuminating  energy  wavelength  divided 
by  the  real  antenna  length  L)  and  R is  the  range  distance  from  the 
antenna  to  the  resolution  cell  on  the  ground.  From  these  two  concepts 
of  resolution  the  number  of  "independent  samples"  which  are  effectively 
combined  to  produce  the  Instantaneous  average  return  power  (PR)  from 
a resolution  cell  can  be  determined: 


. &R  _ 2Rx 
’ L72  ' jT" 


In  a ful ly- focused  synthetic-aperture  radar,  the  azimuth  resolution 
(SR)  is  L/2,  so  there  is  but  one  sample  of  the  random  process  (NS  = 1) 
used  and,  hence  a speckled  appearance  in  the  radar  image11. 


10Bush,  T.  F.  and  F.  T.  Ulaby,  "Fading  Characteristics  of  Panchromatic 
Radar  Backscatter  from  Selected  Agricultural  Targets,"  IEEF  Trans. 
Geosci . Electron.,  Vol . GE-13,  October  1976,  pp.  149-157. 


11  Moore,  R.  K.  and  G.  C.  Thomann,  "Imaging  Radars  for  Geoscience  Use," 
IEEE  Trans.  Geosci,  Electon..  Vol.  GE-9,  July  1971,  pp.  155-164. 
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Averaging  a larger  number  of  "independent  samples"  as  occurs  typically 

in  real-aperture  systems  reduces  signal  fading,  (i.e.,  reduces  the  vari- 

12  11 

ance  and  thus,  smooths  the  image  appearance  ’).  But,  increasing 
the  number  of  samples  averaged  also  degrades  the  final  image  resolution. 

After  detection,  it  can  be  shown  that  the  return  power  from  each 
resolution  cell  is  estimated  by14: 


where  P^  is  the  expected  value  of  the  return  power  from  a resolution  cell 
Y is  a random  variable  with  a standard  chi-square  distribution  having  2NS 
degrees  of  freedom,  and  NS  is  the  number  of  "independent  samples"  averaged 
(Equation  1).  When  the  number  of  independent  samples  being  averaged  is 
large,  (2)  becomes: 


P 


R 


(1  ) 

v/TTs 


(3) 


where  RNis  a normalized  Gaussian  random  variable  with  zero  mean  and  unit 
variance  . 


Zelenka,  J.  S.,  "Comparison  of  Continuous  and  Discrete  Mixed-Integrator 
Processors,"  J.  Opt.  Soc.  Amer. , Vol.  66,  No.  11,  November  1976. 

1 3 

Porcello,  J.  L.,  Norman  G.  Massey,  Richard  B.  Innes,  and  James  M.  Marks, 
"Speckle  Reduction  in  Synthetic-Aperture  Radars,"  J.  Opt.  Soc.  Am., 

Vol.  66,  No.  11,  November  1976. 

14 

Frost,  V.  S.,  J.  L.  Abbott,  V.  H.  Kaupp,  and  J.  C.  Holtzman,  "Derivation 
of  the  Radar  Image  Fading  Characteristics,"  University  of  Kansas 
Center  for  Research,  Inc.,  Remote  Sensing  Laboratory,  Technical  Report 
319-29,  Lawrence,  Kansas,  September  1977. 

^Holtzman,  J.  C. , V.  H.  Kaupp,  J.  L.  Abbott,  V.  S.  Frost,  E.  E.  Komp,  and 
E.  C.  Davison,  "Radar  Image  Simulation:  Validation  of  the  Point  Scatter 
ing  Method,"  Vol.  I,  ETL-0117,  The  University  of  Kansas  Center  for 
Research,  Inc.,  September,  1977,  AD-A053  253. 
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For  terrain  that  can  be  modeled  as  a collection  of  homogeneous,  dis- 
tributed targets  the  value  of  PR  calculated  from  either  (2)  or  (3)  is  the 
best  estimate  of  return  power  after  detection  that  can  be  made  when  each 
resolution  cell  is  treated  as  a single  point.  Even  for  collections  of 
scatterers  not  homogeneously  distributed  this  may  still  be  a good  approach, 
but  statistics  other  than  Rayleigh  may  be  in  order,  for  example,  when 
several  dominant  scatterers  are  located  within  a resolution  cell.  If  the 
terrain  cannot  be  modeled  in  this  way  the  location  and  reradiation  pro- 
perties of  each  of  the  numerous  scattering  points  within  a resolution 
cell  are  modeled,  the  amplitude  and  phase  of  the  return  from  each  point 
is  calculated,  and  the  resulting  phasor  sum  (magnitude)  at  the  antenna  is 
chosen  to  become  the  estimate  of  return  power.  For  high-resolution  and 
synthetic  aperture  radar  applications  of  radar  image  simulation,  either 
approach  is  possible  but  the  latter  is  computationally  inefficient  and 
requires  investment  of  a vast  amount  of  time  to  model  the  ground  properly. 
For  medium  resolution  applications  such  as  those  reported  here.  Equations 
(2)  and  (3)  represent  very  good  estimates  of  the  post-detection  return 
power. 

Usually  after  detection  the  intensity  of  the  video  signal  exiting 
the  receiver  is  recorded  on  film.  This  can  be  expressed  by1: 


D = Ylog10PR  + y 1 og ] qM  ♦ log1QK  (4) 

where  PR  is  given  by  either  Equation  (2)  or  (3);  D is  the  density  of  metallic 
silver  grains  in  the  exposed  transparency  corresponding  to  the  intensity 
of  the  illumination;  K is  a constant  depending  upon  the  exposure  time 
and  the  film  processing  and  development  time;  y is  a positive  constant 
representing  the  slope  of  the  linear  portion  of  the  film  curve  of  density 


Holtzman,  J.  C.,  V.  H.  Kaupp,  J.  L.  Abbott,  V.  S.  Frost,  E.  E.  Komp,  and 
E.  C.  Davison,  "Radar  Image  Simulation:  Validation  of  the  Point  Scatter- 
ing Method,"  Vol.  I,  ETl-0117,  The  University  of  Kansas  Center  for 
Research,  Inc.,  September,  1977,  AD-A053  253. 


versus  logarithm  of  exposure  (the  Hurter-Dri ffield  curve  15  ) where  log^K 
is  the  extrapolated  intercept  of  this  line;  M is  the  transfer  function 
of  the  radar  receiver  (including  all  appropriate  linear  and  non-linear 
effects  such  as  AGC,  Automatic  Gain  Control,  or  saturation)  such  that 
the  video  intensity  Incident  on  the  film  during  exposure  is  specified 
by  MPr  . 

The  next  step  in  the  development  of  the  imaging  model  is  to  define 
how  the  simulation  process  is  to  be  quantized  for  implementing  the  model 
on  a digital  computer.  The  return  power  calculated  for  each  resolution 
cell  is  coded  into  one  pixel  in  the  simulated  image.  Each  pixel  in  the 
image  will  represent  one  precise  shade  of  grey  between  black  (no  return 
power)  and  white  (saturated  signal,  high  return  power).  Thus,  Equation 
(4)  for  the  film  density  must  be  altered  to  reflect  digital  data  process- 
ing requirements.  The  shade  of  grey  (greytone)  for  each  pixel  in  an 
image  can  be  shown  to  be  specified  by  GR  = /,n \ n where  D is  given 


by  Equation  (4),  2n  is  the  number  of  discrete  levels  of  grey  available 
In  a computer  word  having  n bits,  and  g is  the  base  10  logarithm  of  the 
dynamic  range  of  the  radar  signal  being  mapped  into  the  linear  portion 
of  the  film  dyamic  range1.  Using  this  result,  the  final  "greytone" 
equation,  or  imaging  model,  is  given  by:* 


GR  = 


2-1 


y1o9io!e.  + y1°9io!+  lo9lol+Ylog10  <1  + pJ 

P Mc  Kc  V NS 


^ ) \ (5) 
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Goodman,  J.  W. , Introduction  to  Fourier  Optics.  New  York:  McGraw-Hill, 
1968,  pp.  150-153. 

1 Holtzman,  J.  C.,  V.  H.  Kaupp,  J.  L.  Abbott,  V.  S.  Frost,  E.  E.  Komp,.  and 
E.  C.  Davison,  "Radar  Image  Simulation:  Validation  of  the  Point  Scatter- 
ing Method, "Vol.  I,  ETL-0117,  The  University  of  Kansas  Center  for 
Research,  Inc.,  September,  1976,  AD-A053  253. 

*This  final  form  for  the  greytone  equation  is  predicated  upon  using  film 
as  the  final  image  medium.  If  film  is  not  the  final  image  medium, 
the  greytone  equation  would  take  an  appropriate,  alternate  form. 


where  substitution  for  D has  occurred  and  the  result  has  been  calibrated. 
Calibration  is  shown  in  two  phases.  The  first  phase  of  calibration  is 
calculating  the  photographic  density  of  a point  in  the  image  by  the 
expression  2^1  ))o9io[(Ic)t 

signal  exposing  the  film  for  the  calibration  point  (I  = PD  M ) , and  M , 

c k c c 

K , and  y are  the  receiver  transfer  function  and  film  constants  for  the 

calibration  point  as  previously  defined.  The  last  phase  is  setting  G0  , 

n 

the  greytone  to  that  point,  to  a desired  level  (between  0 and  2-1). 
Together  these  calibration  parameters  have  the  effect  of  determining  how 
much  of  the  dynamic  range  of  the  radar  signal  will  be  mapped  into  the 
17-20  dB  dynamic  range  available  in  ordinary  film,  and  exactly  what  portion 
of  the  dynamic  range  of  the  radar  signal  will  be  displayed. 

In  the  development  leading  to  Equation  (5),  the  power  reradiated  as 
a result  of  the  i 1 1 umi nati ng-enerqy/terrain  interaction  was  not  specified. 
An  excellent  model  for  the  return  power  reradiated  from  distributed  targets 
is  given  by  the  radar  equation 

D r2  2 o« 
p . ' o ^ 

R ' -"TT 

(4,)3R  (6) 

where  the  average  transmitted  power  is  represented  by  P,;  the  two-way  gain 
of  the  transmittinq/receiving  antenna  (a  function  of  the  elevation,  and 
azimuth  angles)  is  given  by  G ; the  transmitted  wavelenqth  is  given  by 
the  reflectivity  model,  (a  function  of  wavelenqth  and  local  angle  of  inci- 
dence, among  others)  is  o° ; the  element  of  area  on  the  ground  being  sensed 
(a  function  of  the  ground  slope,  pulse  length,  azimuth  resolution,  and 
altitude)  is  A;  and  the  ranqe  from  the  antenna  to  the  element  of  area  being 
sensed  is  R. 

7Sko1nik,  M.  I.,  (Editor),  Radar  Handbook.  New  York:  McGraw-Hill,  1970. 

g 

Moore,  R.  K.,  "Microwave  Remote  Sensors,"  in  R.  G.  Reeves,  Remote  Sensing 
Manual,  Falls  Church,  Virginia:  American  Society  of  Photogrammetry , 

1975,  Chapter  9. 
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where  I is  the  intensity  of  the 


Certain  conditions  must  exist  for  this  form  of  the  radar  equation  to 
be  valid.  First,  the  area  being  sensed  must  be  a distributed  target.  A 
distributed  target  is  a homogeneous  region  of  a specific  radar  reflectivity 
category  in  which  there  must  be  a large  number  of  individual  scattering 
centers  located  within  a resolution  cell  and  they  must  be  positioned  ran- 
domly 7 . Second,  it  must  be  assumed  that  all  the  parameters  of  the  radar 
equation  are  constant  across  a resolution  cell.  When  these  conditions 
are  satisfied  by  both  the  terrain  being  simulated  and  radar  being  modeled, 
then  Equation  (6)  a particularly  tractable  form  of  the  radar  equation,  can 
be  used  in  the  greytone  Equation  (5)  in  conjunction  with  either  Equation 
(2)  or  (3)  to  estimate  the  return  power  from  each  resolution  cell,  as  was 
done  for  the  results  reported  here. 

1.2.3.  Radar  Geometrical/Propagation  Phenomena 

The  greytone  Equation  (5)  is  the  final  computational  algorithm  of  the 
Point  Scattering  Method  but  it  does  not  represent  the  complete  model.  It 
represents  conversion  of  the  signal  returned  from  each  resolution  element 
into  the  appropriate  grey  shade  for  each  image  pixel  after  the  elevation 
profile,  dielectric  categories,  and  spatial  relationships  of  the  various 
cells  have  been  properly  considered.  The  degree  of  accuracy  of  treatment 
of  these  and  all  the  normal  radar  effects  (such  as  illuminating-energy/ 
terrain  interaction,  layover,  shadow,  range  compression,  etc.)  depends 
upon  the  application  for  which  the  simulated  radar  images  are  being  pre- 
pared. Greater  accuracy  is  required  for  some  applications  than  for  others. 
It  would  be  inefficient  and  expensive  to  treat  all  applications  as  requir- 
ing the  same  detail.  Thus,  while  the  model  is  theoretically  capable  of 
any  degree  of  accuracy  of  treatment,  each  specific  implementation  of  the 
model  is  adjusted  to  meet  the  requirements  of  each  application. 


7Skolnik,  M.  I.,  (Editor),  Radar  Handbook.  New  York:  McGraw-Hill,  1970. 
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The  condition  called  shadow,  manifested  as  a black  area  on  an  image, 
is  caused  by  obstructions  in  the  terrain  which  prevent  illumination  of 
certain  areas  and,  thus,  power  will  not  be  reradiated  continuously  there- 
by causing  gaps  in  the  returned  signal.  Layover  is  also  a result  of 
geometrical/propagation  phenomena  whereby  the  return  power  from  an  object 
nearer  to  the  radar,  by  virtue  of  being  taller  than  its  surroundings,  is 
recorded  sooner,  and  therefore,  closer  to  the  near-range  edge  of  the  image 
than  the  return  from  surrounding  portions  of  the  terrain.  Local  angle  of 
incidence  is  the  effective  angle  made  between  a line  drawn  from  the  anten- 
na to  the  element  of  area  being  sensed  and  the  normal  to  the  element  of 
area  (resolution  cell)  thereby  accounting  for  the  local  slope  (or  tilt)  of 
the  resolution  cell.  This  angle  is  necessary  for  calculating  the  value  of 
o°  for  each  resolution  cell  (Equation  (6)).  If  the  radar  system  being 
modeled  records  in  a slant-range  mode,  then  the  effect  of  range  compres- 
sion must  be  included.  Range  compression  is  the  condition  where  there  is 
a monotonic  decrease  (in  the  range  direction)  of  the  scale  toward  the 
near-range  edge  of  the  image.  If  the  radar  system  being  modeled  records 
in  a ground-range  mode  the  terrain  is  assumed  to  be  flat  and  the  scale  of 
the  image  is  altered  by  the  inverse  of  the  slant-range  scale  compression 
to  preserve  geometric  fidelity.  The  ground-range  mode  adequately  preserves 
geometric  fidelity  for  areas  in  which  relief  is  small  compared  to  the  scale 
of  the  image,  however,  it  accentuates  layover  and  positional  errors  with 
these  becoming  intolerable  in  areas  of  high  relief. 

Treatment  of  all  of  these  effects,  and  more,  is  included  in  the  soft- 
ware implementation  of  the  PSM.  Development  of  the  specific  algorithm 
implemented  for  each  of  the  geometrical /propagation  effects  is  beyond  the 

scope  of  this  report.  Full  development  of  the  simulation  model  is  reported 
1 2 

by  Holtzman,  et  al.  ’ The  total  software  package  computes  the  greytone 
on  a resolution  cell  basis  for  each  pixel  in  the  image  and  accounts  for  all 
geometric/propagation  effects.  There  is  no  one-to-one  correspondance  between 

1 Holtzman,  J.  C.,  V.  H.  -Kaupp,  J.  L.  Abbott,  V.  S.  Frost,  E.  E.  Komp,  and 
E.  C.  Davison,  "Radar  Imaqe  Simulation:  Validation  of  the  Point  Scatter- 
ing Method,"  Vol.  I,  ETL-0117,  The  University  of  Kansas  Center  for 
Research,  Inc.,  September,  1977,  AD-A053  253. 

^Holtzman,  J.  C.,  V.  H.  Kaupp,  J.  L.  Abbott,  V.  S.  Frost,  E.  E.  Komp,  and 
C.  Davison,  "Radar  Image  Simulation:  Validation  of  the  Point 

Scattering  Method," Vol.  II,  ETL-0118,  The  University  of  Kansas  Center 
for  Research,  Inc.,  Sept.  1977,  AD-A053  240. 
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elements  in  the  data  base  and  the  pixels  in  the  simulated  image  since  the 
algorithms  developed  to  model  geometric/propagation  effects  perform  the 
same  function  as  the  radar  system  being  simulated. 


f 


: 

\ 
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1.2.4  Ground  Model 


A model  of  the  ground  must  be  constructed  (ground  truth  data  base) 
for  each  site  for  which  simulated  radar  images  are  to  be  produced.  The 
data  base  is  a symbolic  representation  of  the  different  radar  reflectivi- 
ty categories  and  the  elevation  variations  of  the  terrain  present  in  a 
target  scene.  Since  the  PSM  has  been  implemented  on  the  digital 
computer,  the  data  base  must  be  in  a digital  format. 

Typically  , the  data  base  consists  of  a digital  matrix  containing  four 
dimensions.  These  four  dimensions  are  the  range  and  azimuth  coordinates, 
elevation,  and  radar  reflectivity  (backscatter)  category  of  each  point 
in  the  scene.  It  is  this  matrix  together  with  the  backscatter  data 
(reflectivity  model)  upon  which  the  simulation  program  operates  to  calcu- 
late such  parameters  as  return  power  level,  look-direction,  range,  angle 
of  incidence,  shadow,  layover,  range  compression,  etc. 

The  data  base  can  be  considered  to  be  a symbolic  digital  model  of 
the  physical  (geometric)  and  radar  reflectivity  (dielectric)  properties 
of  the  ground.  Accurate  construction  of  this  digital  terrain  model  is 
crucial  to  the  overall  simulation  effort.  Ground  truth  data  bases  are 
built,  typically,  using  manual  cartographic  feature  extraction  techniques. 
Once  the  source  imagery  (photos,  maps,  etc.)  have  been  obtained  for  a 
particular  site,  a radar/photo-interpreter  uses  these  input  sources  of 
intelligence  data,  and  his  knowledge  and  intuition  to  construct  a back- 
scatter category  map  of  the  area. 

When  this  symbolic  data  base  map  has  been  finished,  it  is  a symbolic 
line  drawing  of  the  boundaries  separating  distributed  targets.  For 
use  on  the  computer,  this  line  drawing  must  be  digitized  and  formed  into 
a matrix.  Digital  elevation  data  must  then  be  added  to  the  category 
matrix  to  produce  the  final  ground  truth  data  base.  After  the  category 
and  elevation  data  are  merged  into  a digital  matrix,  the  ground  truth 
data  base  is  ready  for  input  to  the  radar  simulation  program.  Construction 
details  for  the  data  bases  used  to  produce  the  results  presented  later 
in  this  paper  are  reported  by  Holtzman,  et  al.  (1977)  [2]. 


2 

Holtzman,  J.  C. , V.  H.  Kaupp,  J.  L.  Abbott,  V.  S.  Frost,  E.  E.  Komp,  and 
E.  C.  Davison,  "Radar  Image  Simulation:  Validation  of  the  Point  Scattering 
Method,"  Vol.  II,  ETL-0118,  The  University  of  Kansas  Center  for  Research, 
Inc.,  September,  1977,  AD-A053  240. 


The  choice  of  spacing  of  the  elements  of  the  ground  truth  data  base 
is  governed  by  three  factors.  First,  a terrain  Nyquist  sampling  interval 
for  elevation  data  can  be  determined  by  the  maximum  rate  of  change  of 
elevation  in  the  scene  of  interest.  Second,  each  backscatter  category 
data  point  must  represent  a homogeneous  region  containing  a large  number 
of  randomly  located  scatterers  for  validity  of  the  radar  equation  (6). 
However,  since  the  ultimate  resolution  of  a fully-focused  radar  is 
L/2  (L  is  the  real  antenna  length)  we  need  not  have  the  sample  interval 
<<  L/2.  Thus,  there  are  several  criteria  which  set  lower  bounds  on 
spacing  (i.e.,  minimum  spacing)  between  sample  points  of  the  terrain  in 
the  ground  truth  data  matrix.  On  the  other  hand,  the  maximum  desirable 
spacing  of  sample  points  is  set  by  the  application. 

1.2.5  Reflectivity  Model 

For  each  reflectivity  category  included  in  the  ground  model  (data 
base)  the  dielectric  properties  must  be  modeled.  The  PSM  uses  both 
empirical  backscatter  data  (differential  scattering  cross 
section,  o°)  and  theoretical  results  to  do  this.  The  backscatter  data 
(c°)  model  the  illuminating  energy  and  terrain  interaction  and  specify 
the  fraction  of  energy  reradiated  from  the  ground  back  to  the  antenna. 
Categories  are  represented  in  the  data  base  as,  for  example.  Category  1, 
or  2,  or  3,  or  . . . The  simulation  program,  upon  reading  Category 
1,2,3,  etc.,  for  a point,  solves  the  complex  geometry  relating  the  radar 
platform  to  the  orientation  of  the  ground  at  that  point,  selects  the 
appropriate  backscatter  data-set  for  that  category  and  then  solves  the 
o°  versus  angle-of-incidence  relationship  to  calculate  the  return  power 
via  equation  (5).  Naturally,  the  return  power  calculations  must  reflect 
the  real  situation  in  which  atmospheric  scattering  and  attenuation  may 
become  significant.  In  this  way  the  greytone  is  computed  for  each  pixel 
in  an  image. 


Empirical  a0  versus  angle-of-incidence  data  are  used  in  the  PSM 
wherever  possible  for  the  categories  of  scatterers  typically  encountered 
in  the  terrestrial  envelope  of  interest  for  radar  simulation.  These 
empirical  o°  data  are  taken  both  from  an  extensive  agricultural/soil 
moisture  data  bank  under  development  at  the  Remote  Sensing  Laboratory 
as  well  as  from  the  literature10,16. 


Bush,  T.  F.  and  F.  T.  Ulaby,  "Fading  Characteristics  of  Panchromatic  Radar 
Backscatter  from  Selected  Agricultural  Targets,"  IEEE  Trans.  Geosci. 
Electron.,  Vol . GE- 13,  October  1976,  pp.  149-157. 


Cosgriff,  R.  L.,  W.  H.  Peake,  and  R.  C.  Taylor,  "Terrain  Scattering  Proper- 
ties for  Sensor  System  Design,"  Eng.  Experiment  Station  Bulletin, 

181,  Vol.  29,  Ohio  State  University,  Columbus,  Ohio,  May,  1960. 


The  feasibility  of  using  the  PSM  for  making  reference  radar  images 
for  a terminal  guidance  concept  using  the  Correlatron  was  investigated 
in  a joint  effort  between  the  RSL  (Remote  Sensing  Laboratory,  The  Univer- 
sity of  Kansas)  and  the  ETL  (U.S.  Army  Engineer  Topographic  Laboratories, 
Fort  Bel  voir,  Virginia).  Figure  2 is  reproduced  from  reference  (1)  to 
illustrate  the  apportionment  of  work  between  RSL  and  ETL  because  the 
same  plan  was  followed  in  this  study.  The  three  major  classifications 
of  work  at  RSL  are:  (1)  produce  data  bases,  (2)  obtain  terrain  reflectivi- 
ty data,  and  (3)  produce  reference  scenes.  The  activities  at  ETL  included 
supplying  digital  elevation  data,  producing  reference  scene  film,  and 
running  Correlatron  tests. 

The  first  and  most  time  consuming  task  accomplished  was  the  construc- 
tion of  the  ground  truth  data  base.  The  Pickwick  Landing  Dam  site  had 
been  previously  selected  for  the  lower  altitude  scene  produced  earlier 
on  the  basis  of:  (1)  existence  of  real  radar  video  data,  (2)  diversity 
of  radar  scattering  types  in  the  scene,  (3)  availability  of  source  data 
to  construct  the  ground  truth  data  base,  (4)  accessibility  of  empirical 
backscatter  coefficient  data  to  describe  the  scatterers  within  the  site. 
Likewise,  the  Pickwick  site  met  these  criteria  for  the  higher  altitude 
scene  reported  in  this  work.  Two  data  bases  were  constructed  for  the 
reference  radar  images  formed.  The  first  data  base  was  constructed  to 
support  reference  scenes  for  the  two  larger  scale  (lower  altitude)  images 
(called  bands  1 and  2)  and  one  for  the  two  smaller  scale  (higher  altitude) 
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images  (called  bands  3 and  4).  Part  of  the  work  of  making  the  data  base 
for  the  two  larger  scale  images  was  done  and  reported  earlier  (a  data  base 
sufficient  for  the  band  1 image  was  produced  and  reported  in  reference  (1)). 
Construction  of  these  data  bases  is  discussed  in  a subsequent  section. 

Terrain  reflectivity  data  were  sought  for  the  microwave  reflectivity 
categories  built  into  the  two  Pickwick  ground  truth  data  bases.  As  a conse- 
quence of  the  poorer  spatial  resolution  built  into  the  band  3 and  4 data  base, 
quite  often  many  distinctly  different  scattering  regions  were  grouped 
together.  However,  since  the  planimetry  source  data  (aerial  photographs) 
displayed  much  finer  resolution,  it  was  possible  to  determine  the  constitu- 
ency of  areas  where  a mixture  of  targets  would  be  grouped  as  one  composite 
scattering  type.  There  are  no  clear  cut  methods  on  how  to  weigh  the  importance 
of  the  response  of  individual  scattering  regions  (of  sub-resolution  cell 
size)  within  a resolution  cell.  Therefore,  it  was  a venture  into  the 
unknown,  guided  by  intuition,  that  led  to  the  combination  of  empirical 
backscatter  data  that  were  utilized  for  the  bands  3 and  4 data  base  cate- 
gories for  inhomogeneous  resolution  cells.  Weighting  the  response  of  scat- 
tering regions  by  the  relative  percentage  of  area  coverage  within  a resolution 
cell  was  often  used. 

Upon  completion  of  these  activities  (construction  of  the  data  bases 
and  cataloguing  appropriate  backscatter  data  from  diverse  sources),  the 
three  reference  scenes  were  generated  and  stored  on  computer-compatible 
magnetic  tape  for  shipment  to  ETL. 

The  completed  reference  scenes  were  sent  to  ETL  where  the  reference 
scene  film  was  produced  from  the  digital  tapes.  The  film  containing 
the  reference  scenes  was  stored  on  ETL's  Correlatron  test  stand  and  was 
tested  versus  "live  " radar  video  data  that  had  been  recorded  via  a video 
recorder  in  an  earlier  flight-test  program. 
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2.0  REFERENCE  SCENE  GROUND  TRUTH  DATA  BASES 

Two  ground  truth  data  bases  have  been  constructed  to  support  making 
reference  radar  scenes  for  the  Pickwick  Landing  Dam  target  site.  One 
data  base  was  prepared  for  the  bands  1 and  2 scenes  (largest  scale  and 
lowest  altitude  scenes)  and  one  for  the  bands  3 and  4 scenes  (smallest 
scale  and  highest  altitude  scenes).  These  data  bases  were  developed 
to  model  the  terrain  in  the  topoqraphic  region  of  the  target  site  in  the 
states  of  Alabama,  Tennessee,  and  Mississippi  centered  on  the  northwest 
corner  of  the  power  house  at  the  dam  (coordinates  34“  04'  15"  N by 
88°  15'  05"  W). 

Each  of  the  data  bases  is  a digital  model  of  the  terrain  (digital 
matrix)  in  the  target  site  containing  the  range  and  azimuthal  coordinates 
of  each  point,  the  microwave  backscatter  cateqory  of  each  point,  and  the 
elevation  above  sea-level  of  each  point.  The  bands  1 and  2 data  base 
is  discussed  in  section  (2.1)  and  the  bands  3 and  4 data  base  in  section 
(2.2). 

2 . 1 Bands  1 and  2 Data  B ase 

Construction  techniques,  source  data,  and  details  for  the  bands  1 
and  2 data  base  are  adequately  related  in  reference  (17).  The  discussion 
will  not  be  repeated  here,  however,  specification  of  pertinent  parameters 
regarding  this  data  base  is  presented  in  Table  1. 

Table  2 lists  the  microwave  backscatter  categories  determined  to  be 
in  the  Pickwick  target  site  within  the  area  comprising  the  bands  1 and  2 
data  base.  The  actual  backscattering  (o')  data  for  each  category  used  in 
conjunction  with  this  data  base  are  listed  in  Appendix  A. 
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Table  1:  Parameter  Specifications:  Bands  1 and  2 Data  Base  (Pickwick  Site). 


Data  Base  Size:  Square  of  approximately  19.8  km  (12.3  miles)  on  a side; 

Site  Location:  Center  located  at  coordinates  34°  04'  1 5"N  by  88°  15'  05"W 
(Pickwick  Landing  Dam,  Tenn.); 

Spatial  Sample  Size:  6.25m  (20.5  feet)  in  both  range  and  azimuth; 

Elevation  Data  Accuracy:  Estimated  to  be  approximately  + 3m  ( + 10  feet); 

Backscatter  Category  Resolution:  Estimated  to  be  approximately  30.5  m 

(100  feet)  in  both  range  and  azimuth; 

Number  of  Elements  in  Digital  Matrix:  10,042,561  (3169  samples  at  6.25  m 

per  sample  in  both  range  and  azimuth) 
where  N-3169  records  and  M=  3169  points; 

Scale:  1:24,000; 

Source  Intelligence  Data  Used: 

Elevation  Data:  Provided  by  ETL  and  produced  by  DMA 
(Defense  Mapping  Agency)  from  20  foot  contour  data  of 
USGS  (United  States  Geological  Survey)  standard  7 1/2' 
quadrangle  maps; 

Category  Data:  Spatial  geometry  and  detail  was  obtained 

from  USGS  standard  7 1/2'  quadrangle  maps  (1:24,000  scale) 
and  distributed  category  boundaries  interpreted  from 
1:100,000  high-resolution  aerial  photographs  (Note: 
the  geometry  of  the  terrain  in  a standard  7 1/2'  USGS 
quadrangle  is  displayed  as  a modified  polyconic). 
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Table  2 


Pickwick  Test  Site  Microwave 

Category  Identification 
Number 

Backscatter  Categories:  Bands  1 and  2 

Data  Base 

Category  Description 

100 

Group  Targets  - Roads 

110 

Heavy  duty  improved  roads  (none  present) 

120 

Medium  duty  improved  roads 

130 

Light  duty  improved  roads 

140 

Unimproved  roads 

200 

Group  targets  - Railroads 

230 

Fish  Pond  Dikes 

240 

Water  Plant  Plumbing 

300 

Group  Targets  - Water  Bodies 

310 

Small  Impoundments 

320 

Small  streams  and  rivers 

350 

Large  streams  and  rivers 

360 

Large  Impoundments  (Pickwick  Lake) 

400 

Group  targets  - Forested  areas 

450 

Wooded  Marshes 

500 

Group  Targets  - Forested  Areas 

600 

Group  Targets  - Agricultural  Lands 

610 

Bare  ground 

620 

Soybeans 

630 

Corn 

640 

Milo 

650 

Wheat 

660 

Orchards 

670 

Garden  plots 

700 

Group  Targets  - Grass-covered  areas 

(parks,  etc. ) 

800 

Pickwick  Dam 

810 

Blockhouse 

820 

Spi 1 lway 

900 

Small  buildings 

910 

Large  buildings 

Note:  The  microwave  reflectivity  categories  listed  in  this  table  represent 
a level  of  spatial  detail  of  approximately  18  meters. 
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2.2  Bands  3 and  4 Data  Base 


Specification  of  pertinent  parameters  regarding  the  data  base  is 
presented  in  Table  3.  A listing  of  the  microwave  backscatter  categories 
determined  to  be  in  the  Pickwick  target  site  within  the  bands  3 and  4 
data  base  area  is  arrayed  in  Table  4 and  the  actual  backscatter  (0°)  for 
each  category  used  in  conjunction  with  the  data  base  are  included  in 
Appendix  B. 

The  bands  3 and  4 data  base  was  constructed  by  following  the  basic 
philosoohy  developed  for  the  bands  1 and  2 data  base.  The  only  signifi- 
cant differences  in  the  construction  of  the  bands  3 and  4 data  base  from 
that  reported  earlier  for  bands  1 and  2^ 7 are  the  source  data  used. 

The  bands  1 and  2 data  base  was  developed  from  standard  7 1/2' 

USGS  (United  States  Geological  Survey)  Quadrangle  maps  at  a scale  of 

1:24,000  (the  scene  geometry  and  elevation  data  were  obtained  from  such 

maps  and  the  backscatter  region  bouridaries  were  interpreted  from  1:100,000 

scale  aerial  photographs).  The  data  base  for  bands  3 and  4 was  developed 

1 8 

from  an  orthophoto  which  itself  had  been  constructed  via  the  UNAMACE 
computer  system.  The  orthophoto  is  a controlled  mosaic  of  high-resolution 
aerial  photographs  in  which  the  geometry  is  rectified  to  become  a tangent 
plane  approximation  to  the  Earth  at  a specified  point.  This  signifies 
a basic  difference  between  the  bands  1 and  2 and  bands  3 and  4 data  bases: 
The  geometry  of  the  former  is  a modified  polyconic  projection  according 
to  standard  USGS  mapping  practices  while  the  geometry  of  the  latter  is  a 
tangent  plane  projection  centered  on  the  specified  Pickwick  target.  In 
the  bands  1 and  2 data  base  both  the  geometry  of  the  scenes  and  elevation 
data  for  the  scene  were  obtained  from  a modified  polyconic  projection 
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Table  3 Parameter  Specifications:  Bands  3 and  4 Data  Base 

(Pickwick  Site) 


Data  Base  Size:  Square  of  approximately  74.5  km  (46.3  miles)  on  a side; 

Site  Location:  Center  located  at  coordinates  34n  04'  15"  N by  88°  15'  05"  W 
(Pickwick  Landing  Dam,  Tenn.); 

Spatial  Sample  Size:  25  m (82.0  feet)  in  both  range  and  azimuth; 

Elevation  Data  Accuracy:  Estimated  to  be  anproximately  + 10m  (+32.8  feet); 

Backscatter  Category  Resolution:  Estimated  to  be  approximately  244  m (800  feet) 

in  both  ranqe  and  azimuth; 

Number  of  Elements  in  Digital  Matrix:  8,880,400  (2,980  samples  at  25  m per 

sample  in  both  range  and  azimuth)  where 
N=2980  records  and  m=3169  points. 

Source  Intelligence  Data  Used: 

Elevation  Data:  Provided  by  ETL  from  the  output  of  the  UMAMACE 
elevation  data  computer  program; 

Category  Data:  Spatial  neometry  and  detail  was  obtained  from 

1:100,000  scale  orthophoto  (rectified  geometry) 
and  distributed  category  boundaries  were  interpreted 
from  the  high-resolution,  1:100,000  scale  aerial 
photographs  which  also  served  as  the  source  for  the 
orthophoto.  (Note:  The  geometry  of  the  terrain  in 
the  orthophoto  was  rectified  for  a tangent  plane 
approximation  to  the  Earth  centered  on  the  target 
center) . 
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whereas  in  the  bands  3 and  4 data  bases  both  the  geometry  and  elevation 
data  were  obtained  from  a tangent  plane  projection. 

A simple  reason  explains  this  difference  in  projections  between  the 
large  scale  (bands  1 and  2)  and  small  scale  (bands  3 and  4)  data  bases: 
Expediency.  When  the  original  work  was  performed  resulting  in  the  band  1 
data  base,  it  was  decided  to  utilize  the  USGS  maps  in  the  modified  polyconic 
projection  instead  of  the  tangent  plane  projection  of  the  UNAMACE  ortho- 
photo because  the  geometry  contained  in  the  maps  was  believed  to  be  a better 
approximation  to  reality  at  the  target  site  than  the  tangent  plane  offered 
t by  UNAMACE.  In  addition,  a group  of  researchers  at  ETL  were  using  the  UNAMACE 

data  as  source  information  and  the  use  of  USGS  maps  provided  an  evaluation 
of  an  alternate  source.  However,  when  the  bands  3 and  4 work  was  started 
the  UNAMACE  data  (elevation  and  category)  was  selected  as  source  material 
because  elevation  data  were  available  from  ETL  at  a desirable  sampling  rate 
(one  sample  every  25m)  and  a 1:100,000  scale  orthophoto  was  available,  whereas 
the  elevation  data  corresponding  to  the  USGS  maps  were  arrayed  at  too  fine 
a sampling  rate  (6.25  m/cell)  across  more  than  fifty  (50)  separate  computer- 
compatible  magnetic  tapes,  each  tape  having  a unique  coordinate  system.  In 
view  of  the  fact  that  in  an  operational  environment  source  data  such  as 
high-resolution  aerial  photographs  would  be  used  for  construction  of  reference 
scenes,  it  was  decided  to  use  the  UNAMACE  data  (which  uses  aerial  photographs 
as  source  material)  for  the  smaller  scale  data  base  because  of  its  availability 
and  to  save  rectifying  the  coordinate  systems  of  fifty  tapes  in  order  to 
select  desired  elevation  data  from  them. 

The  sequence  followed  when  constructing  a data  base  (either  bands  1 and 
2 or  bands  3 and  4)  consists  of  several  basic  steps: 

1)  Obtain  elevation  data  at  a suitable  scale; 

2)  Obtain  source  data  for  the  category  map; 

3)  Develop  the  category  map; 

a)  Interpret  the  source  data  and  draw  the  various  features 
and  boundaries  onto  a stable-base  drawing  medium; 

b)  Digitize  the  category  maps; 

c)  Develop  matrix  from  the  digitized  form  of  the  maps 
which  has  a category  specified  for  each  of  the  elements; 

UK 
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4)  Merge  the  elevation  data  and  category  matrix  into  the 
digital  matrix  representing  the  final  data  base; 

5)  Rotate  the  data  base  into  the  required  orientation. 

2.2.1  Elevation  Data  for  Bands  3 and  4 Data  Case:  Pickwick  Site 

1 ft 

Elevation  data  prepared  via  the  UNAMACE  software  system  were 
obtained  from  ETL.  These  data  were  stored  on  a single  computer-compatible 
magnetic  tape.  A sample  of  the  elevation  was  stored  for  each  25  m cell 
in  the  tangent  plane  projection  used  by  the  UNAMACE  system.  The  accuracy 
of  these  data  were  estimated  to  be  approximately  + 10  m. 
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2.2.2  Source  Data  for  Bands  3 and  4 Data  Base:  Pickwick  Site 


Source  data  for  the  Pickwick  site  consisting  of  high-resolution 
aerial  photographs  (1:100,000  scale,  nominally)  and  an  orthophoto 
(1:100,000  scale  tangent  plane  projection  of  target  site  prepared  by 
the  UNAMACE  system)  were  obtained  from  ETL.  Only  these  source  data  were 
used  for  the  bands  3 and  4 data  base;  no  radar  or  other  ancillary  data 
were  used. 


2.2.3  Develop  the  Bands  3 and  4 Category  Map:  Pickwick  Site 

The  bands  3 and  4 data  base  was  constructed  at  a scale  of  approxi- 
mately 1:150,000.  This  was  accomplished  via  a scale-reduction-table 
and  was  done  to  accomodate  requirements  of  the  digitizer  system  that  was 
to  be  used.  The  1:100,000  scale  orthophoto  was  placed  on  the  table 
and  a projection  of  the  scene  was  focused  upon  a glass  plate,  forming 
a backlit  image  on  the  stable-base  material  which  had  been  placed  over 
the  glassplate.  The  planimetry  of  the  data  base  was  traced  at  a scale 
of  1:150,000  directly  onto  the  drawing  material  from  this  backlit  image. 

The  planimetry  of  the  data  base  was  developed  via  standard  feature 
extraction  technique.  A photo  interpreter  defined  the  boundaries  out- 
lining homogeneous  regions  (homogeneous  at  radar  wavelengths)  and  drew 
these  on  the  stable-base  drafting  material.  Care  was  exercised  through- 
out construction  of  this  feature  map  to  hold  spatial  resolution  to  approx- 
imately 0.05  inches  on  the  orthophoto  which,  at  1:150,000  scale,  corresponds 
to  approximately  800  feet  on  the  qround.  This  means  that  any  radar  image 
constructed  from  this  data  base  will  not  contain  scene  variations  of  less 
than  800  feet;  a scattering  region  cannot  be  di fferent iated  from  its 
surroundings  unless  it  has  a spatial  extent  of  800  feet,  or  more.  This 
is  the  inherent  resolution  designed  into  the  data  base  and  is  not  necessarily 
the  resolution  of  any  resultant  rada--  reference  scene. 

The  planimetry  of  the  Pickwick  site  was  drawn  on  two  sheets  of  the 
stable-base  drawing  material:  One  of  the  sheets  contained  the  homogene- 
ous regions  of  two-dimensional  extent  (areas  having  a finite  width  and 
length)  and  the  other  sheet  contained  the  "category  map",  or  "distributed 
target  nap".  The  first  of  these  was  called  the  "distributed  category  map" 
and  the  second  was  called  the  "cultural  target  map",  or  "hard  target  map". 

On  the  "category  map"  was  drawn  the  boundaries  of  each  different  back- 
scattering  region  such  as  water  bodies,  forests,  agricultural  fields,  etc., 
identified  to  be  within  the  tarqet  site.  On  the  "cultural  target  map"  was 
drawn  all  the  linear  features  such  as  roads,  railroads,  etc.  The  locations 
of  buildings  were  marked  on  the  "cultural  target  map"  as  was  one  type  of 
area  target:  Cities. 

4-  i 
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Table  4 Pickwick  Test  Site  Microwave  Backscatter  Categories:  Bands  3 & 4 

Data  Base 


Category  Identification* 
Number 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 


Category  Description 
Forested  Area 

Agricultural  Lands  - Undifferentiated 

Reservoir 

River 

Low  Grass  Area 

Agricultural  Lands  - Flooded 

Roads  - Heavy  Duty 

Roads  - Medium  Duty 

Roads  - Light  Duty 

City  Category  - Mostly  Structures 

" - Mostly  Trees 

" - Mixed  Structures 

& Trees 

" - Open  Land 

Utility  Cut  - Low  Grass 


! 


* Category  1 was  not  used. 


Cities  were  modeled  as  consisting  of  different  regions,  each  region 
having  the  average  backscattering  properties  of  one  of  four  (4)  city 
backscatter  categories.  The  four  (4)  city  backscatter  categories  and 
the  o°  values  used  for  the  production  of  radar  reference  images  are  listed 
in  Table  5.  As  can  be  seen  from  the  table  this  listing  of  four  (4)  city 
backscatter  categories  is  very  subjective.  Nothing  quantitative  is  being 
claimed  about  the  radar  echo  from  cities  or  their  backscatter  properties. 
The  goal  of  this  work  was  to  develop  a data  base  for  producing  reference 
radar  scenes  for  terminal  guidance.  Reference  scenes  must  be  as  omni- 
directional as  possible  because  no  a priori  knowledge  exists  concerning 
the  direction  of  approach  of  the  projectile.  This  being  the  case,  it 
was  desired  to  model  the  cities  in  the  target  region  in  a manner  which 
would  acknowledge  their  presence  without  separately  accounting  for  each 
building,  road,  intersection,  corner,  etc.  (this  information  would  be 
desirable  for  a simulation  having  highly  directional  properties).  It 
was  decided  to  do  more  than  just  "paint  the  city  regions  white,"  thus 
the  four  (4)  separate,  subjective  city  backscatter  categories  were  created 
for  the  data  base  in  an  attempt  to  refine  the  city  into  a "reasonable" 
number  of  distinct  levels. 

After  both  the  "category  map"  and  the  "cultural  target  map"  were 
finished  (all  of  the  appropriate  regions  and  1 i near  features  of  the  Pick- 
wick site  were  drawn  on  one  of  the  two  sheets  of  drafting  material), 
the  line  content  of  the  drawings  was  digitized.  A large  table  digitizer 
interfaced  to  a minicomputer  was  used  for  this  task.  Each  sheet  was  separ- 
ately placed  on  the  digitizer  table  and  an  operator  first  specified  a 
coordinate  system  (the  coordinate  system  specified  was  that  of  the  UNAMACE 
elevation  data)  and  then  traced  the  data  using  the  tables  cursor.  As  the 
operator  traced  a line,  the  system  periodically  sampled  the  location  of 
the  cursor  with  respect  to  a wire  grid  built  into  the  table.  As  the 
operator  encountered  a new  region,  or  line,  he  identified  it  to  the  mini- 
computer by  typing  its  label  via  a teletypewriter.  The  minicomputer  then 
converted  each  sample  position  into  an  x,y  pair  of  numbers  in  the  specified 
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Table  5 City  Backscatter  Categories 


Category 

No. 

Description  of  Category 

Backscatter  Values 

1 

Mostly  hard  targets  such  as  a 
business  district  in  a city; 

+ 20  dB  isotropic* 

2 

Primarily  hard  targets  but  also 
including  a significant  propor- 
tion of  grassy  or  forested 
regions; 

- 5 dB  isotropic* 

3 

Primarily  grassy,  or  forested 
region  but  also  including  a 
significant  proportion  of  hard 
targets ; 

- 9 dB  isotropic* 

4 

Mostly  grassy,  or  forested 
regions  in  a city  such  as  a 
park. 

Low  Grass (Category 
16  in  Appendix  B) 

* The  backscatter  values  cited  are  nominal  values  attained  in  the  final 
simulations.  To  attain  these  values,  a cosine  rule  was  utilized  to 
modify  the  data  for  input  to  the  reference  scene  simulation  program. 
Copies  of  plots  of  these  data  are  included  in  Appendix  B. 
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coordinate  system  and  stored  these  data  with  the  label  for  ultimate  output 
on  computer-compatible  magnetic  tape.  This  portion  of  the  data  base  con- 
struction task  was  inordinately  long  because  of  numerous  problems  encounter- 
ed such  as  hardware  failures  (the  cursor  became  defective  causing  a length y 
delay  while  the  problem  was  investigated  andrepai red) , schedule  problems 
(both  the  system  was  overscheduled  as  well  as  understaffed),  operator 
training  problems  (a  new  operator  was  hired  and  trained  by  digitizing  this 
data  base-  a deplorable  state  of  affairs  but  one  over  which  we  had  no 
control  as  the  system  on  which  this  work  was  being  done  was  the  only  one 
in  the  local  area),  etc.  The  final  product  of  this  digitizing  effort  was 
computer-compatible  magnetic  tapes  on  which  were  stored  the  data  from  both 
maps,  separately. 

After  both  maps  were  digitized,  these  digital  boundary  data  were  con- 
verted into  a matrix  which  has  a category  specified  for  each  of  its  ele- 
ments. In  principle  this  task  is  not  difficult  because  it  is  only  desired 
to  sort  the  boundary  data  by  their  x and  y values  and  "fill-in"  the  correct 
category  for  each  element  in  the  matrix  bounded  by  a unique  line.  In 
practice  this  is  not  as  trivial  as  it  might  sound  because  the  data  base 
matrix  contained  ten  million  (10^)  elements,  there  being  several  hundred 
thousand  boundary  values,  alone,  and  it  was  highly  desirable  to  program 
the  computer  to  do  this  task  without  external  manipulations  of  the  data. 

The  programs  developed  to  accomplish  this  task  have  been  discussed  pre- 
viously,^ so  they  will  not  be  discussed  here. 

This  task  was  complicated  by  errors  created  during  the  digitizing 
of  the  maps  by,  in  particular,  the  lack  of  experience  on  the  part  of  the 
table  operator.  These  errors  caused  an  extensive  delay  in  producing  the 
final  category  matrix  because  they  required  human  investigations  to  deter- 
mine the  nature  of  the  problems  and  to  identify  and  to  implement  corrections. 
It  would  be  highly  desirable  to  create  the  finished  category  matrix  while 
the  digitizing  is  being  done,  and  to  provide  a visual  reconstruction  of 


1 9 

Holtzman,  J.  C.,  V.  H.  Kaupp,  J.  L.  Abbott,  V.  S.  Frost,  E.  E.  Komp,  and 
E.  C.  Davison,  "Radar  Image  Simulation:  Validation  of  the  Point  Scatter- 
ing Method,"  Vol.  I,  ETL-0117,  The  University  of  Kansas  Center  for  Re- 
search, Inc.,  September,  1977,  AD-A053  253,  pp.  131-163. 


37 


each  step  of  the  process  via  a graphics  terminal  (or  equivalent)  so  that 
errors  could  be  dispensed  with  while  the  map  is  still  on  the  table  and  while 
the  digitizer  operator  can  assess  what  he  views  as  the  finished  product. 

After  both  the  "category  map"  and  the  "cultural  target  map"  were 
developed  into  completely  specified  matrices,  they  were  merged  into  one 
array.  The  data  from  the  "cultural  target  map"  was  given  priority  over 
the  data  from  the  "category  map"  thus  creating  a single  matrix  containing 
both  area  targets  and  linear,  or  cultural  (buildings,  etc.)  targets. 

2.2.4  Merge  Elevation  and  Category  Data:  Bands  3 and  4 

As  the  category  data  had  been  digitized  in  the  same  coordinate  system, 
with  the  same  cell  size,  as  the  UNAMACE's  elevation  data,  it  was  easy 
to  merge  the  two  (elevation  and  category  data).  One  complication  arose 
in  this  task,  the  UNAMACE  elevation  data  were  anomalous  over  very  flat 
regions  such  as  water  surfaces.  This  was  solved  by  specifying  the  eleva- 
tion for  water  bodies;  the  elevation  of  the  river  was  set  at  362  feet  and 
of  the  reservoir  at  414  feet  above  mean  sea  level.  This  correction  obvious- 
ly is  inadequate  wherever  the  river  or  reservoir  blends  smoothly  into  a 
flat  flood  plain,  for  there  the  anomalous  elevation  data  still  exists. 

The  percentage  of  the  target  site  for  which  this  condition  exists  was  so 
small  the  problem  was  not  corrected. 

After  completion  of  the  task,  the  bands  3 and  4 data  base  was  finished; 
it  was  ready  for  the  production  of  radar  reference  scenes. 

2.2.5  Rotate  Bands  3 and  4 Data  Base  into  Desired  Orientation 

The  radar  reference  scenes  which  were  to  be  formed  from  this  data 
base  must  have  a specified  orientation  with  respect  to  true  north.  The 
first  pixel  stored  on  the  computer-compatible  magnetic  tape  of  the  reference 
scenes  arrayed  in  a rectangular  grid  format  must  be  the  northwest  corner 
of  the  scene.  The  UNAMACE  coordinate  system  was  skewed  30°  from  true 
north.  Rather  than  rotate  the  data  base  at  this  time,  it  was  decided 
to  rotate  after  the  simulations  (reference  radar  images)  were  produced 
and  while  the  data  were  still  in  a polar  coordinates  system.  This  tech- 
nique was  considerably  easier  to  implement  and  was  much  less  expensive  of 
computer  time. 
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3.0  RADAR  REFERENCE  SCENE  SIMULATION  MODEL 


No  changes  in  the  Point  Scattering  Method  (PSM)  implementation 

for  producing  radar  reference  scenes  were  required  for  bands  2,  3,  and  4 

reference  scenes.  The  software  implementation  developed  for  the  band  1 

20 

reference  scene  and  previously  reported  was  adequate  for  the  bands 
2,3,  and  4 scenes.  The  implementation  previously  developed  is  summar- 
ized here  for  convenience. 

The  PSM  was  to  be  adapted  to  model  the  terminal  guidance  problem 
consisting  of  a PPI  (Plan-Position  Indicator)  radar  and  Correlatron. 

The  first  step  in  this  special izati on  of  the  PSM  to  a specific  system 
was  to  attempt  to  describe  the  operating  parameters  of  the  PPI  radar, 
itself.  However,  limited  information  about  the  operating  characteris- 
tics of  both  the  radar  and  Correlatron  was  available.  Therefore,  in  the 
absence  of  system  design  data,  the  reference  scene  simulation  software 
was  developed  assuming  an  ideal  system.  For  instance,  the  PPI  radar 
(for  simulation  purposes)  was  given  constant  azimuthal  gain  between  its 
3 dB  points  and  no  sidelobes  (an  aspiration  for  any  antenna  designer.'). 
The  elevation  pattern  was  chosen  to  be  (cs c:'tf)  (cos h),  where  ts  is  the 
depression  angle.  Past  the  rf  portion,  the  receiver  of  the  ideal  system 
was  made  to  map  the  received  power  into  video  intensity.  A realistic 
film  transfer  characteristic  was  employed  (logarithmic)  with  a linear 
dynamic  range  of  20  dB.  Outside  this  range,  either  in  the  "toe"  or 
"shoulder"  of  the  exposure  curve,  lack  of  sufficient  exposure  or  satura- 
tion, respectively,  would  result. 

It  was  secondly  considered  whether  there  should  be  modifications  to 
the  reference  scene  simulation  model  to  account  for  the  Correlatron. 
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It  was  assumed  to  have  identical  paths  for  both  the  reference  scene  and 
real  video.  The  process  of  converting  a reference  scene  stored  on  photo- 
graphic film  to  a video  signal  was  assumed  to  be  linear.  And,  of  course, 
identical  tests  run  at  different  times  were  assumed  to  result  in  the  same 
degree  of  cross  correlation.  All  of  these  criteria  were  assumed  to  be  the 
case  for  the  Correlatron. 

Complicating  the  situation  of  specializing  the  PSM  for  constructing 
reference  scenes  was  the  fact  that  the  direction  of  approach  of  the 
vehicle  was  not  specified.  To  optimize  the  chances  of  high  correlation, 
and  to  allow  the  reference  scenes  to  be  useful  despite  the  radar  position 
and  angle  of  approach,  it  was  necessary  to  make  the  reference  scenes  as 
nearly  omnidirectional  as  possible.  This  was  shown  to  dictate  a nadir- 
looking antenna  because  of  the  angular  dependence  of  both  radar  shadow 
and  the  backscattered  fields.  The  only  information  about  the  system 
available  before  either  constructing  the  data  bases  or  the  reference 
scenes  was:  (1)  the  reference  scene  altitude  and  (2)  the  corresponding 
diameter  of  each  simulated  PPI  image.  Thus,  each  image  was  formed  with 
the  radar  centered  over  the  power  house  of  the  Pickwick  Landing  Dam,  look- 
ing radially  outward,  as  though  its  trajectory  was,  at  least  momentarily, 
vertical  to  the  earth. 

Figure  3 illustrates  the  peculiar  image  format  of  the  guidance  PPI 
radar  being  employed,  in  comparison  with  the  ordinary  PPI  radar  scan 
format.  Data  are  recorded  by  the  test  radar  for  a full  circular  sweep 
of  the  scene  instead  of  the  usual  sector  associated  with  PPI  radars.  The 
terrain  imaged  by  the  radar  beam  is  within  an  annular  ring  bounded  at  the 
near  range  by  35°  (incidence  angle)  and  at  the  far  range  by  65°<  The 
reference  scene  simulation  model  does  not  produce  imagery  in  exactly  this 
format  because  of  (1)  the  likelihood  of  centering  and  angle  of  approach 
errors,  and  (2)  the  use  of  the  Correlatron  as  the  diagnostic  device. 

The  direction  of  approach  of  the  projectile  and  real  PPI  radar  and 
the  center  of  the  real  imagery,  being  unknown  before  forming  reference 
scenes,  might  occur  any  place  within  a circular  region  of  the  target 
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Figure  3.  Special  PPI  Image  Format. 


depending  upon  ballistic  guidance  variables.  This  presented  a formidable 
problem  for  the  reference  scene  software  because  the  reference  scene  had 
to  be  capable  of  a high  degree  of  cross  correlation  with  the  real  image 
centered  anywhere  in  the  reference  scene  region  so  that  guidance  corrections 
could  be  developed.  According  to  information  received  from  ETL,  the 
angle  of  the  tangent  to  the  trajectory  (as  measured  from  the  vertical) 
at  each  specified  altitude  would  be  at  most  45°.  To  ascertain  that  the 
"live"  video  data  would  fall  well  within  the  reference  scene  at  the  max- 
imum angle  of  approach  (measured  from  vertical)  an  allowance  was  made 
which  would  enlarge  the  reference  scene  radius.  These  considerations 
together  dictated  that  the  simulated  image  boundary  extend  to  about  75°. 
These  conditions  (uncertainty  of  approach  and  center  of  the  real  PPI 
video  signal)  imposed  necessary  conditions  on  the  simulation  effort: 

(1)  No  "holes"  were  allowed  in  the  reference  scene  even  though  one 
existed  in  the  real  data;  (2)  The  reference  scene  must  be  larger  than 
the  real  PPI  scene  to  accommodate  "centering"  errors;  (3)  Angles  of  inci- 
dence of  the  real  data  were  in  the  range  35°  to  65°,  meaning  all  of  the 
reference  scene  should  also  be  in  that  range  even  though  the  actual  geo- 
metry of  the  reference  scene  would  decree  a range  from  0°  to  approximately 
75°. 

The  difficulties  caused  by  look-direction  errors  could  be  of  severe 
magnitude.  If  the  PPI  radar  approaches  the  target  (northwest  corner  of 
the  power  house)  off  course,  then  look-direction  errors  between  real  and 
and  reference  scenes  between  0°  and  130°  occur.  Look  direction  effects 
are  most  significant  in  ground  scenes  having  considerable  local  relief 
because  the  direction  and  length  of  shadows  (and  layover)  in  radar  images 
are  determined  by  the  look  direction.  Figure  4 illustrates  the  problem 
for  a look  direction  error  of  180°  between  the  real  and  simulated  radar 
image.  Test  sites  having  significant  local  relief  variation  would  appear 
very  different  depending  upon  the  direction  of  approach  of  the  real  PPI 
radar.  Fortunately,  the  Pickwick  test  site  had  only  a modest  amount  of 
local  relief,  so  this  is  one  problem  which  was  not  tackled.  But  it  certain- 
ly warrants  attention  if  a test  site  having  a considerable  amount  of  local 


Figure  4.  Comparison  of  Reference  Scene  to  PPI  Radar  Image  Format. 


relief  is  selected  at  some  point  in  the  future. 

Other  aspects  of  the  direction  problem  were  treated  and  their  impacts 
minimized  in  the  software  implementation  of  the  simulation  model.  As 
should  be  obvious  from  Figure  4,  as  the  look  direction  changes  so  also 
does  the  incidence  angle  change.  As  is  well  known,  ground  radar  return 
data  (o°)  for  the  same  target  varies  by  many  decibels  over  the  range 
0°  to  70°  incidence  angle.  For  this  reason,  the  reference  scene  simulation 
software  could  not  be  set  up  to  produce  a simulated  image  according  to 
the  actual  geometry  of  the  problem.  If  this  were  done,  even  if  the  65° 
circle  happened  to  fall  always  on  the  same  category  and  thus  would  be 
a constant  shade  of  grey  in  the  real  image,  the  same  65°  circle  would  trace 
out  a path  on  the  reference  image  which  could  conceivably  vary  from 
black  to  white. 

In  the  range  of  incidence  angles  in  the  real  image,  35°  to  65°,  most 
radar  ground  return  curves  are  relatively  linear  and  have  relatively 
shallow  slopes  (nothing  factual  or  quantitative  implied  here,  this  is 
a qualitative  arguement).  The  antenna  function  in  the  range  direction 
(look  direction)  over  this  sane  range  of  incidence  angles  tends  to  compen- 
sate for  the  slope  of  the  o°  data,  producing  for  a number  of  ground  return 
categories  a relatively  uniform  return,  thereby  minimizing  the  problem 
caused  by  different  look  directions.  For  these  reasons,  it  was  decided 
that  the  minimum  angle  of  incidence  in  the  simulated  reference  scene 
would  be  35°  and  the  maximum  would  be  65°.  The  area  in  the  reference 
scene  lying  between  35°  and  65°  angle  of  incidence  was  simulated  normally. 
The  area  lying  within  the  35°  circle  was  imaged  as  though  the  angle  of 
incidence  was  a constant  35°.  And,  the  area  lying  outside  the  65°  circle 
was  all  simulated  as  if  it  were  65°.  This  solution  did  not  attempt  to 
model  the  real  video  exactly,  but  rather  did  attempt  to  minimize  discrep- 
ancies between  the  reference  scene  and  the  "live"  video  produced  in  flight. 
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This  is  not  to  jay  that  local  slope  variations  were  not  accounted  for; 
they  were  indeed,  incorporated.  What  is  meant  is  that  the  incidence  angle 
(e)  between  the  antenna  "boresight"  and  the  local  vertical  was  always 
in  the  range  35°  to  65°.  Local  slope  variations  then  altered  the  incidence 
angle  to  the  local  incidence  angle  (e^),  just  like  in  the  simulation 
models  described  in  References  [1]  and  [2].  In  fact,  the  limitations 
imposed  on  minimum  and  maximum  values  of  come  strictly  from  the 
local  relief  in  the  scene. 

This  solution  to  the  angle  of  incidence  problem  created  data  handling 
problems  for  the  computer  program,  and  data  base  problems.  For  instance, 

35°  angle  of  incidence  specifies  a resolution  cell  size  for  short  pulse 
and  narrow  beamwidth  radars.  Yet,  the  geometry  of  the  data  base  indicates 
that  as  data  base  cells  get  closer  to  the  center  (in  polar  coordinates), 
they  get  larger  in  the  range  direction  and  smaller  in  the  azimuth  direction. 
This  problem  was  minimized  by  accurately  modeling  another  feature  of  the 
real  PPI;  it  recorded  data  in  ground  range  mode.  Ground  range  mode  means 
that  (for  a flat  earth)  equal  size  objects  located  in  the  near  and  far 
range  will  have  equal  sizes  in  the  image  format.  This  is  normally  ac- 
complished by  applying  a nonlinear  sweep  to  the  electron  beam  of  the 
viewing  CRT.  But  for  simulation  purposes,  it  simply  meant  building 
the  simulation  data  base  with  equal  size  cells  in  the  range  direction. 

It  should  be  noted  at  this  time  that  in  the  presence  of  terrain  having 
significant  relief,  ground  range  mode  introduces  large  distortions,  a 
fact  to  keep  in  mind  for  such  future  sites. 

In  summary,  the  general  point  scattering  radar  image  simulation  model 
was  specifically  tailored  to  the  special  requirements  imposed  to  simulate 
reference  scenes  for  use  on  the  Correlatron.  The  software  implementation 
of  the  reference  scene  simulation  model  included  the  following  special 
features. 
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Table  6 Reference  Scene  PSM  Special  Features 


(1)  360°  PPI  image  scan  format 

(2)  Simulated  area  was  larger  than  the  real  image  to  allow 
"centering"  errors. 

(3)  No  holes  allowed,  the  reference  scene  was  completely  filled-in 
with  radar  image  simulations. 

(4)  Minimum  angle  of  incidence  = 35° 

(5)  Maximum  angle  of  incidence  = 65° 

(6)  Local  angle  of  incidence  was  properly  treated. 

(7)  In  the  reference  scene,  the  area  between  0°  and  35°  was  simulated 
at  a constant  35°  angle  of  incidence. 

(8)  The  area  between  35°  and  65°  was  simulated  normally. 

(9)  The  area  outside  65°  was  simulated  at  a constant  65°. 

(10)  Variations  due  to  angle  of  incidence  difference  between  real 
and  simulated  image  were  minimized. 

(11)  Reference  scenes  were  formed  in  the  ground  range  model. 

(12)  Layover  and  shadow  were  properly  included. 

(13)  Local  slope  variations  in  the  terrain  were  properly  included. 
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4.0  TERRAIN  BACKSCATTER  DATA:  PICKWICK  SITE 
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The  Point  Scattering  Method  (PSM)  for  radar  image  simulation  utilizes 
empirical  backscatter  (o°)  data  (whenever  possible),  in  combination  with 
relief  data,  to  model  the  radar  echo  from  terrain.  The  backscatter  data 
(o°)  supplies  the  needed  information  about  the  type  of  terrain  and  pro- 
vides the  mechanism  to  predict  what  tone  will  be  recorded  for  that  type 
of  terrain  in  an  image.  The  elevation  data  determines  the  radar-to-target 
distance  and  the  direction  of  the  local  normal  to  the  surface  and,  thus, 
determines  the  angle  between  the  radar  antenna  and  the  normal  to  the 
surface.  Together,  the  o°  and  elevation  data  are  used  to  predict  the 
tone  for  each  pixel  (picture  element)  in  a simulated  radar  image;  the 
o°  data  predict  what  percentage  of  transmitted  energy  illuminating  each 
piece  of  ground  will  be  reradiated  back  to  the  radar  receiving  antenna 
for  the  particular  range  and  local  angle  of  incidence  considerations 
that  exist. 

After  specification  of  the  PPI  radar  configuration  (polarization  and 
frequency)  and  after  the  different  types  of  scatterers  had  been  recognized 
for  the  Pickwick  site,  efforts  were  begun  to  find  appropriate  terrain 
return  data.  Targets  were  first  grossly  classified  as  distributed  or 
cultural  and  symbolic  representation  was  necessary  for  the  latter. 

Distributed  targets  are  homogeneous  regions  with  the  same  microwave 
scattering  properties  throughout  the  extent  of  a resolution  cell.  Each 
homogeneous  region  must  be  at  least  as  large  as  the  resolution  element 
of  the  radar  being  modeled,  the  individual  scattering  centers  must  be 
randomly  located,  and  there  must  be  a large  number  of  scattering  centers 
in  each  resolution  cell  within  a homogeneous  region.  When  these  conditions 
are  satisfied,  an  average  value  of  the  scattering  cross-section  (a°)  can 
be  used  to  model  the  radar  return  from  these  homogeneous  areas  of  terrain 
(distributed  targets).  Most  of  the  terrain  located  in  the  reference 
scene  data  base  of  the  Pickwick  site  satisfied  these  criteria.  Thus, 
differential  scattering  cross-section  data  (0°)  were  used  to  model  the 
radar  return  properties  of  the  terrain  in  the  reference  scenes  formed  from 
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the  Pickwick  data  base.  Actually,  empirical  a°data  were  used  (as  opposed  | 

to  theoretical).  These  data  were  obtained  from  the  literature  and  from 

the  RSL  data  bank.  The  best  match  that  could  be  found  between  available 

empirical  a°data  and  the  identified  distributed  targets  in  the  Pickwick 

site  was  sought  (the  data  used  are  presented  in  Appendices  A and  B). 

Cultural  targets  are  here  defined  to  be  manmade  objects  and  features. 

Their  radar  returns  are  characterized  by  specular  reflection.  They  cannot 
be  modeled  as  distributed  targets.  That  is,  the  return  from  hard  targets 
are  not  readily  predicted  by  the  greytone  equation.  To  model  cultural 
targets  by  digital  computer  is  an  exceptionally  complex  task  requiring 
tremendous  detail  about  each  such  target  to  be  included  in  the  data  base. 

An  evaluation  was  made  of  the  kinds  of  cultural  targets  present  in  the 
Pickwick  site,  their  orientations,  and  the  number  of  them.  It  was  con- 
cluded that  symbolic  modeling  of  these  targets  would  provide  the  best 
probability  of  correlation.  By  symbolic  modeling  we  mean  that  the  location 
of  each  cultural  target  was  pin-pointed  in  the  data  base  but  the  orientation, 
size,  geometry,  etc.,  was  not.  Cultural  targets  were  assumed,  for  the 
purposes  of  reference  scene  formation,  to  be  isotropic  radiators  having 
a constant  effective  differential  scattering  cross-section.  This  most 
certainly  is  not  accurate.  But,  consider  the  task:  Cross  correlation 
with  "live"  video  data  with  an  unknown  heading  and  center  with  respect 
to  the  center  of  the  reference  scene.  The  orientation  (not  the  location) 
of  cultural  targets  was  unknown.  Given  no  a priori  knowledge  of  the 
heading  of  the  "live"  data,  it  would  not  be  possible  to  properly  simulate 
all  the  corner  reflectors,  etc.,  of  cultural  targets  to  match  the  look- 
direction  of  the  "live"  data.  The  best  that  could  be  done  would  be  to 
accurately  mark  the  location  of  them  and  assume  them  to  be  isotropic 
radiators.  Then,  given  a direction  of  approach  for  the  real  data, 
those  cultural  features  properly  aligned  would  be  in  the  real  data  which, 
of  course,  would  match  very  nicely  with  the  simulation.  Those  not  properly 
aligned  would  be  some  lower  shade  of  gray  which  should  still  improve  cor- 
relation. Thus,  it  was  concluded  that  symbolic  simulation  of  cultural 
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features  would  enhance  cross  correlation  between  the  real  and  simulated 
images,  and  this  is  the  way  they  were  treated  when  forming  reference 
scenes. 

Since  two  (2)  different  data  bases  were  constructed,  one  for  bands 
1 and  2,  and  one  for  bands  3 and  4,  the  definition  of  backscatter  targets 
identified  to  be  in  the  sites  was  specialized  to  take  advantage  of  the 
different  resolutions.  The  backscatter  data  used  for  each  data  base 
are  discussed  in  the  following. 

4. 1 Backscatter  Data  for  Bands  1 and  2:  Pickwick  Site 

The  backscatter  data  used  to  make  bands  1 and  2 reference  radar 
scenes  are  presented  in  Appendix  A.  The  data  are  identified  as  to 
category,  frequency,  polarization,  and  source  from  which  the  data  were 
acquired. 

4.2  Backscatter  Data  for  Bands  3 and  4:  Pickwick  Site 

The  backscatter  data  used  to  make  bands  3 and  4 reference  scenes, 
except  the  city  backscatter  data  which  are  listed  in  Table  5,  are 
presented  in  Appendix  B.  The  data  are  identified  as  to  category, 
frequency,  polarization,  and  source  from  which  the  data  were  acquired. 


5.0  RADAR  REFERENCE  SCENE  FORMATION 

It  is  to  this  point  where  all  the  previous  work  has  been  aimed: 
Production  of  radar  reference  scenes.  The  various  kinds  of  input  data 
have  been  assembled  in  the  above  work  and  they  are  ready  to  be  combined 
appropriately  into  reference  scenes.  Figure  5 illustrates,  in  block 
diagram  form,  how  the  separate  pieces  fit  together  and  the  interactions 
that  occur  between  them  to  make  reference  scenes. 

Figure  5 shows  the  flow  of  data,  from  start  to  finish,  when  forming 
a radar  reference  scene.  The  left  side  of  Figure  5 illustrates  where 
the  data  comes  from,  and  the  right  side  where  it  goes.  The  central 
portion  of  Figure  5 shows  the  formation  of  reference  scenes. 

As  can  be  seen  from  Figure  5,  the  radar  reference  scene  implementa- 
tion of  the  PSM  is  for  a digital  computer.  Three  separate  computer 
programs  (separate  according  to  function)  must  be  run  sequentially  in 
order  to  form  radar  reference  scenes.  The  first  of  these  computer 
programs  (POLAR  CONVERSION)  converts  a rectangular  ground  truth  data  base 
matrix  into  a polar  radar  data  base  array  as  dictated  by  the  desire  to 
simulate  a polar-scanning  radar  (PPI)  for  making  reference  scenes.  The 
second  computer  program  (REFERENCE  SCENE)  accepts  the  polar- 'orm  radar 
data  base  from  the  first,  solves  all  the  complicated  geometrical  relation- 
ships between  the  radar  (located  at  a specific  point)  and  each  resolution 
element  on  the  ground,  computes  the  power  returned  from  each  resolution 
element  to  the  radar,  incorporates  the  appropriate  antenna  correlation 
function  between  resolution  elements,  and  produces  an  output  array  of 
image  density  values  (called  greytones).  The  third  program  (RECTANGULAR 
CONVERSION)  converts  the  array  of  greytones  which  were  output  from  the 
second  program  in  a polar  array  into  a rectangular  grid  matrix  having 
either  a format  compatible  with  viewing  the  finished  simulation  for 
evaluation  purposes  or  a format  compatible  with  the  Correlatron  test 
equipment  utilized  at  ETL.  Communication  between  programs  is  accomplished 
via  computer-compatible  magnetic  tapes.  These  tapes  contain  intermediate 
products  and  only  serve  as  temporary  storage  either  for  error  analysis  or 
for  input  to  the  next  program,  thus  their  formats  will  not  be  discussed  here. 


In  the  followinq  five  (5)  sections,  radar  reference  scene  formation 
is  summarized.  First,  the  work  preliminary  to  making  radar  reference 
scenes  is  briefly  recapitulated  (section  5.1).  Second,  third,  and 
fourth,  computer  programs  used  for  forming  radar  reference  scenes  are 
discussed  (sections  5.2,  5.3,  and  5.4,  respectively).  Fifth,  mechanisms 
for  evaluation  of  radar  reference  scenes  are  briefly  mentioned  (section 
5.5). 


5. 1 Preliminary  Work 


As  can  be  seen  from  Figure  5,  the  first  tasks  to  be  accomplished, 
in  the  chain  of  events  which  leads  ultimately  to  radar  simulation  for 
reference  scenes,  are  specification  of  the  guidance  system  to  be  modeled, 
and  the  terrestrial  location  of  the  target  site.  After  system  and  site 
have  been  specified  three  important  tasks  can  be  performed:  1)  An  imple- 
mentation can  be  designed  of  the  PSM  which  is  specialized  for  the  system 
and  its  unique  features;  2)  A ground  truth  data  base  can  be  developed  which 
is  a digital  representation  of  the  terrain  in  the  target  site;  and,  3) 

A catalogue  of  o°  data  can  be  obtained  which  is  a model  for  the  dielectric 
properties  of  the  terrain  backscatter  categories  identified  in  the  target 
site.  Upon  completion  of  these  tasks,  the  data  necessary  to  satisfy  input 
requirements  of  the  computer  proqrams  which  form  the  reference  radar  scenes 
are  avai lable. 

An  implementation  of  the  PSM  was  designed  around  the  navigation  system 
consisting  of  PPI  radar  and  Correlatron.  This  implementation  was  a 
specialization  of  the  PSM  for  making  reference  radar  scenes  for  terminal 
guidance.  The  details  of  this  special ization,  the  work  performed,  and 
the  resultant  reference  scene  model  are  summarized  in  section  (3.0).  The 
resultant  model  is  a set  of  programs,  called  REFERENCE  SCENE  SIMULATION, 
which  has  been  tailored  specifically  to  make  reference  scenes  for  the 
Correlatron  guidance  system. 

Two  different  scale  ground  truth  data  bases  were  developed  for  a target 
site  centered  on  the  Pickwick  Landing  Dam,  Tennessee.  One  data  base  was 
developed  to  support  formation  of  reference  scenes  for  the  two  (2)  large 
scale  (low  altitude)  scenes,  and  the  other  for  the  two  (2)  small  scale 
(high  altitude)  scenes.  The  development  of  these  data  bases  is  summarized 
in  section  (2.0).  The  completed  data  bases  were  rectangular  grid  matrices 
of  N records,  each  record  containing  M points,  each  point  being  eighteen 
(18)  bits  long.  N and  M are  specified  in  Tables  1 and  3 for  the  two  dif- 
ferent data  bases.  The  first  twelve  (12)  bits,  of  each  eighteen  (18)  bit 
point  in  a data  base,  contain  the  elevation  in  feet  above  mean  sea  level 
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for  that  point,  and  the  trail inq  (6)  bits  are  used  to  identify  up  to 
sixty-three  (63)  different  backscatter  categories.  Each  data  base  is 
stored  on  a 9-track,  1600  bpi  computer-compatible  magnetic  tape  in 
a binary  unformulated  format.  These  two  ground  truth  data  bases  stored 
on  magnetic  tape  serve  as  fundamental  inputs  to  the  reference-radar- 
scene-  format  ion  computer  programs. 

Terrain  backscatter  (a°)  data  were  obtained  and  catalogued  for  all 
the  different  reflectivity  targets  in  the  Pickwick  site.  Thirty  (30) 
different  categories  were  identified  for  the  bands  1 and  2 data  base 
(Table  2)  and  fourteen  (14)  for  the  bands  3 and  4 data  base  (Table  4). 
This  work  and  the  data  obtained  are  discussed  in  section  (4.0).  These 
data  were  obtained  in  the  form  of  backscatter  versus  angle  of  incidence 
[ct° ( e ) ] . Third-order  polynomials  were  fit  to  these  data  according  to 
standard  Least-Squares  techniques  and  the  coefficients  were  catalogued 
for  input  to  the  simulation  program.  These  data  are  lead-card  input 
to  the  program  for  ease  of  changing  categories  and  backscatter  data. 

They  are  stored  in  a P x 4 matrix  where  P represents  the  number  of 
categories  and,  thus,  data  sets,  for  a data  base,  each  category  having 
four  (4)  entries,  there  being  one  entry  for  each  coefficient  of  the 
third-order  polynomial.  Validity  of  using  a third-order  polynomial  to 
represent  the  average  properties  of  the  j°  data  can  be  observed  in 
Appendices  A and  B. 

Upon  completion  of  these  tasks,  the  input  data  for  the  reference 
scene  computer  programs  have  been  obtained  (ground  truth  data  base  and 
terrain  backscatter)  and  the  PSM  model  has  been  specialized  to  represent 
the  guidance  system.  The  next  task  is  to  construct  radar  reference 
scenes  from  the  input  data. 


5.2  Polar  Conversion 


The  first  step  performed  by  the  sequence  of  computer  programs 
illustrated  in  Figure  4 is  POLAR  CONVERSION,  the  computer  program 
which  converts  the  rectangular  grid  matrix  of  the  ground  truth  data 
base  Into  a polar  array  which  is  compatible  with  the  scanning  format 
of  the  PPI  radar  being  modelled.  POLAR  CONVERSION  actually  consists 
of  three  (3)  distinct  computer  programs,  POLAR  CREATE,  POLAR  ARRAY, 
and  ARRAY  FIX;  copies  of  these  programs  are  provided  In  Appendices 
C,  D,  and  E,  respectively. 

Three  programs  have  been  developed  to  perform  this  straight-forward 
function  of  rectangular-to-polar  coordinates  conversion  in  order  to 
minimize  the  costs  of  performing  it.  The  first  program,  POLAR  CREATE, 
is  a highly  computational  program  requiring  minimal  core  storage.  This 
program  accepts  as  input  the  ground  truth  data  base  stored  on  digital 
magnetic  tape,  computes  the  polar  address  (r,e)  of  each  point  as  it 
Is  read  from  tape,  performs  one-dimensional  compression  on  the  data  and 
stores  the  data  sequentially  on  an  intermediate  magnetic  tape.  Compres- 
sion arises  from  the  fact  that  the  rectangular  version  (original  version) 
of  the  data  base  contains  finer  resolution  and,  consequently,  more 
sample  points  than  are  desired  in  the  polar  data  base.  The  number  of 
elements  desired  in  the  final  polar  data  base  (calculated  from  the  radar 
resolution  parameters)  Is  a control  parameter  input  to  the  program 
and  is  used  to  quantize  the  polar  conversion  calculations  and,  thus, 
produce  a mechanism  for  compression.  Table  7 lists  the  number  of  ele- 
ments and  the  resolution  element  size  (estimated  to  represent  one 
independent  sample;  NS=  1,  equation  (l))  utilized  for  each  of  the  four 
(4)  bands  of  reference  scenes  produced  here. 


Table  7 


Re  sol  >.  j t i Ele-ent  Size:  Bands  1-4 


Resolution  Cel'  Size* 


Reference 
Scene  Band  Mo. 

r 

Mo.  of  Elements 
in  Polar  Data  Base 

1 

30.5m 

(100 

feet) 

1/2° 

105,120 

2 

30.5m 

(100 

feet) 

1/2° 

210,240 

3 

100  m 

(328 

feet) 

1/2° 

128,160 

1 

100  n 

(328 

feet) 

1/2° 



255,600 

1- 

♦Resolution  in  the  ranpe  direct-,  ! aint.ained  constant  at  the  value 

specified,  as  was  resolution  in  tin  i.-ii  jt  I rection. 


The  second  program,  POLAR  ARRAY,  is  a core- intensive  but  computational ly- 
minimal  program.  This  program  accepts  as  input  the  intermediate  magnetic 
tape  output  from  POLAR  CREATE,  orders  the  data  from  the  sequential  rec- 
tangular file  to  the  correct  polar  array,  and  stores  the  POLAR  ARRAY 
on  an  intermediate  magnetic  tape.  Two-dimensional  compression  is  per- 
formed by  this  proqram.  This  follows  the  same  philosophy  as  does  the 
earlier  compression.  The  category  data  are  compressed  on  a priority 
basis;  the  highest  priority  category  brought  from  the  rectangular  data 
base  for  each  polar  cell  is  retained,  the  others  deleted.  Elevation  data 
are  averaged.  The  elevation  value  stored  in  any  polar  cell  is  the 
average  of  the  elevation  of  all  the  rectangular  points  which  were  mapped 
into  it. 

The  third  program,  ARRAY  FIX,  was  created  to  rectify  the  problem 
created  by  the  fact  that  the  rectangular- to-polar-conversion  mapping  is 
less  than  one-to-one  in  the  center  portion  of  the  polar  data  base.  This 
program  interrogates  the  nearest  filled  neighbor  to  determine  the  category 
and  elevation  of  a polar  cell  found  empty.  Upon  satisfactory  operation 
of  this  program,  the  completed  polar  data  base  is  stored  on  an  intermediate 
magnetic  tape. 

In  this  way  POLAR  CONVERSION  functions  to  convert  a large  rectangular 
grid  matrix  data  base  into  a smaller  polar  array  data  base  for  minimal  cost. 
Cost  is  minimized  because  the  amount  of  core  (a  high-cost  component)  which 
must  be  used  during  the  computational ly-intensive  portions  of  the  operation 
is  reduced  to  the  bare-essential  amount.  The  output  of  POLAR  CONVERSION 
is  a computer-compatible  magnetic  tape  containing  the  ground  truth  data  base 
arrayed  in  a polar  format  with  the  correct  resolution  to  support  directly 
the  REFERENCE  SCENE  programs,  the  next  programs  which  must  be  run. 

{ . 
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5.3  Reference  Scene 


Reference  to  Figure  4 will  show  the  second  step  performed  is 
REFERENCE  SCENE,  the  computer  program  which  actually  performs  the 
simulation  of  the  guidance  radar  system  and  forms  the  desired  simulations 
of  radar  images.  REFERENCE  SCENE  actually  consists  of  two  (2)  computer 
programs,  POWER,  and  GREYTONE:  copies  of  these  programs  are  provided  in 
Appendices  F and  G,  respectively. 

Two  programs  have  been  developed,  instead  of  one,  to  minimize  the 
costs  and  improve  the  operational  efficiency  of  running  the  programs 
many  times.  The  first  program,  POWER,  is  a computational ly-intensi ve 
program  requiring  minimal  core  storage.  This  program  accepts  as  input 
both  the  polar  ground  truth  data  base  on  digital  magnetic  tape  from 
program  POLAR  CONVERSION  and  the  terrain  backscatter  data  lead-card 
input,  calculates  the  average  power  exiting  the  receiver  on  a pixel-by- 
pixel  basis  for  each  pixel  (picture  element,  previously  called  "point") 
in  the  final  scene,  and  stores  these  data  on  an  interim  magnetic  tape. 

POWER  recognizes  each  radial  record  in  the  polar  data  base  as  the  scan 
line  corresponding  to  the  energy  returned  from  one  pulse  of  the  radar, 
each  point  in  the  record  corresponds  to  a resolution  element.  The  polar 
data  only  contains  a record  for  successive,  independent  scan  lines  (pulses). 

The  PRF  (Pulse  Repetition  Frequency)  of  a radar  system  is  normally 
quite  high  with  successive  pulses  producing  a return  having  a large  over- 
lap with  several  preceding  pulses.  POWER  calculates  a new  scan  line  (new 
pulse)  of  data  only  for  scan  lines  which  are  independent  of  one  another 
(they  do  not  overlap  each  other),  and  calculates  the  dependency  of  over- 
lapping pulses  statistically  as  noted  by  equations  (2)  and  (3)  together 
with  their  incorporation  in  the  greytone  equation  (5).  This  is  done  in 
order  to  minimize  both  the  size  of  the  polar  data  base  and  the  computational 
load  (and,  thus,  cost)  required  to  produce  radar  reference  scenes. 

Similarly,  POWER  recognizes  each  point  in  a scan  line  from  the  polar 
data  base  as  an  independent  resolution  element  in  the  radial  direction, 
and  calculates  the  data  relating  each  point  on  the  ground  to  a pixel  in 
the  image.  Dependency  of  overlapping  samples  in  the  radial  direction  is 
statistically  incorporated  in  the  model  (equations  2 and  3). 


For  each  point  In  the  data  base  (each  resolution  element  on  the 
ground),  POWER  solves  the  geometry  relating  the  position  of  the  radar 
platform  (three-dimensional  position)  to  the  point  and  calculates  the 
slope  of  the  terrain,  and  the  angle  of  incidence  and  the  range  between 
platform  and  point.  These  calculations  are  made  sequentially  for  each 
point  of  each  record  as  the  tape  containing  the  polar  ground-truth 
data  base  is  read  into  the  computer.  Upon  determining  these  parameters 
for  a point,  POWER  enters  the  main  computational  algorithm  of  the  program 
which  calculates  the  average  power  exiting  the  receiver  from  that  point. 
This  calculation  of  power  uses  the  slope  of  the  ground  (two-dimensional 
slope),  the  angle  of  incidence  between  radar  and  ground  (both  normal  and 
local  angles),  the  range  from  platform  to  the  point  on  the  ground  being 
interrogated,  the  power  pattern  of  the  antenna,  the  category  identifica- 
tion from  the  polar  data  base,  backscatter  data  from  the  a0  file,  and 
the  transmitter/recei ver/image  model  incorporated  for  the  radar  system 
whose  response  is  being  simulated.  All  of  these  variables  and  parameters 
are  combined  appropriately  for  calculating  the  estimate  of  power  exiting 
the  radar  receiver  for  each  point  on  the  ground.  In  this  way  POWER  cal- 
culates the  average  power  exiting  the  receiver  on  a pixel-by-pixel  basis. 
The  resultant  data  are  stored  on  an  interim  magnetic  tape  for  further 
processing  in  later  stages.  The  data  are  ordered  sequentially  on  this 
tape  in  the  same  form  as  the  polar  array  in  which  the  polar  radar  data 
base  was  input. 

The  second  program  of  REFERENCE  SCENE  was  developed  to  incorporate 
the  spatial  relationships  between  adjacent,  independent  cells  decreed 
by  the  antenna  pattern,  and  to  convert  the  resultant  estimates  of  power 
into  greytones.  This  program,  GREYTONE,  calculates  the  spatial  relation- 
ships between  cells  via  an  autocorrelation.  The  shape  and  length  of 
the  autocorrelation  are  Input  parameters.  Upon  completing  the  auto- 
correlation, GREYTONE  converts  the  data,  power,  into  density  values, 
greytones,  quantizes  them  Into  the  desired  number  of  bits,  and  biases 
the  range  to  that  desired  for  ultimate  storage  in  a photograph. 
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The  bias  required  to  display  a desired  power  range  in  an  image  is  an  input 
parameter  to  GREYTONE.  The  desired  mapping  ratio  of  power  exiting  the 
receiver  into  density  in  the  photograph  is  also  an  input  parameter;  the  por- 
tion of  the  radar  dynamic  range  (XdB)  desired  to  be  mapped  into  the  dy- 
namic range  of  the  photograph  (17-20dB)  is  specified.  Thus,  upon  specifi- 
cation of  an  autocorrelation  function  shape  and  length,  and  quantizing 
parameters  (bias  and  mapping  ratio,  or  gain)  GREYTONE  operates  on  the  power 
map  input  via  digital  magnetic  tape  from  the  previous  program,  POWER, 
producing  the  greytone  map  of  the  final  image  on  a pixel-by-pixel  basis. 

The  greytone  data  are  stored  on  an  intermediate  digital  magnetic  tape. 

The  stored  data  order  is  still  the  same  as  the  input  polar  ground  truth 
data  base. 

At  this  point,  the  radar  image  simulation  work  is  complete  but  the 
data  are  still  stored  in  a polar  array.  The  following  program  converts 
these  data  from  polar  back  to  rectangular  coordinates  for  compatibility 
with  either  standard  raster-scan  format  display  devices  for  evaluation 
or  the  Correlatron  for  testing  purposes. 


5.4  Rectangular  Conversion 
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Reference  to  Figure  4 will  show  the  third  step  performed  is  RECTANGULAR 
CONVERSION,  the  computer  program  which  converts  the  simulated  radar  image 
from  a polar  array  to  a rectangular  grid  matrix.  RECTANGULAR  CONVERSION 
actually  consists  of  two  (2)  computer  programs,  RECTANGULAR  CREATE,  and 
RECTANGULAR  ARRAY;  copies  of  these  programs  are  provided  in  Appendices 
H and  I,  respectively. 

These  two  programs  exist  for  the  same  purpose  as  POLAR  CREATE  and 
POLAR  ARRAY  (section  5.2)  and  perform  the  inverse  operations  of  them. 
RECTANGULAR  CREATE  requires  input  specification  of  the  size  of  the 
rectangular  grid  desired  for  the  output  data.  The  size  of  this  grid  is 
dependent  upon  the  purposes  for  which  the  data  have  been  created.  If 
the  data  have  been  created  for  display  and  evaluation  purposes,  the  size 
of  the  grid  is  entirely  determined  by  the  size  of  the  area  to  be  viewed 
and  the  size  limitations  of  the  display  device.  If  the  data  have  been 
created  for  testing  on  the  Correlatron,  the  size  of  the  output  grid  is 
specified  to  be  921  x 921  pixels.  Upon  specification  of  the  size  of  the 
output  rectangular  grid  matrix,  RECTANGULAR  CREATE  calculates  again  the 
polar  address  (r,e)  of  each  point  in  the  specified  rectangular  array  and 
stores  these  data  on  an  intermediate  magnetic  tape. 

RECTANGULAR  ARRAY  requires  input  specification  of  the  data  format  of 
the  magnetic  tape  to  be  output.  For  display  and  evaluation  purposes,  the 
format  is  defined  to  be  raster  format  with  the  word  length  for  the  grey- 
tones  of  each  pixel  , and  whether  a positive  or  a negative  is  desired,  to 
be  specified.  For  testing  on  the  Correlatron,  the  output  format  is  de- 
fined to  be: 

9 track  digital  magnetic  tape; 

1600  bpl; 

921  records; 

921  pixels  per  record; 

0 Corresponds  to  white; 

255  Corresponds  to  black. 
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Upon  specification  of  these  input  parameters,  RECTANGULAR  ARRAY  converts 
the  data  output  from  RECTANGULAR  CREATE  into  a completely  specified 
rectangular  grid  having  the  desired  output  format,  and  stores  these 
data  on  a digital  magnetic  tape. 

When  these  data  have  been  created  for  testing  on  the  Correlatron 
(the  appropriate  specifications  have  been  made),  the  tape  output 
from  RECTANGULAR  ARRAY  contains  the  desired  reference  radar  scene  for 
one  of  the  bands. 

5. 5 Testing  and  Evaluation 

Upon  completion  of  the  previous  activities  and  having  a data  tape 
output  from  RECTANGULAR  ARRAY,  the  data  are  available  for  testing  and 
evaluation.  Standard  display  devices  such  as  a monochrome  teTevision 
set  which  is  interfaced  to  the  digital  computer  via  IDECS*  are  employed 
at  RSL  for  image  evaluation.  This  capability  is  exploited  frequently 
for  trouble-shooting  problems  and  for  validating  reference  scenes  before 
they  are  sent  to  ETL  for  testing  on  the  Correlatron. 

Once  a simulated  radar  image  nas  been  produced  and  has  passed  the 
evaluation  phase  at  RSL,  the  output  of  RECTANGULAR  CREATE  is  run  again 
through  the  reference  scene  format  version  of  RECTANGULAR  ARRAY  for  pro- 
ducing a digital  tape  having  the  right  format  for  the  Correlatron  test 
system.  This  tape  contains  a reference  radar  image  for  testing  at  ETL 
on  the  Correlatron.  The  tape  is  mailed  to  ETL  and  the  reference  scene 
is  run  through  the  Correlatron  system  by  personnel  at  ETL.  The  results 
of  running  our  reference  radar  scene  against  "live"  video  data  on  the 
Correlatron  are  communicated  to  RSL  from  ETL. 


★ 

IDECS  is  the  acronym  standing  for  Image  Discrimination,  Enhancement, 
Combination,  and  Sampling.  The  IDECS  is  an  analog  image  enhancement 
station  located  at  the  RSL. 
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6.0  PROBLEMS  AND  SOLUTIONS 


6.1  Geometric  Fidel i tv 


When  planning  the  various  aspects  involved  in  the  construction  of 
the  data  base  for  bands  3 and  4,  it  was  decided  to  produce  the  hand-drawn 
planimetry  at  a smaller  scale  than  the  source  intelligence  (1:150,000 
instead  of  1:100,000  scale).  This  was  decided  so  that  the  resultant 
hand-drawn  map  could  be  digitized  on  a digitizer  table  which  was  more 
readily  available  than  the  one  used  for  the  bands  1 and  2 data  base.  To 
accomplish  this  reduction  in  scale,  a reducing  table  was  used.  Use  of 
this  reducing  table  caused  geometric  fidelity  problems. 

The  ortho  photo  which  was  used  as  the  source  intelligence  for  the 
bands  3 and  4 data  base  was  too  large  to  fit  completely  on  the  reducing 
table.  This  meant  that  only  a small  segment  of  the  data  base  could  be 
developed  for  a given  orientation  of  the  ortho  photo  on  the  table.  Then 
both  the  ortho  photo  and  the  hand  drawing  had  to  be  moved  and  re-registered. 
Not  apparent  at  the  time  of  construction  but  discovered  later,  this  moving 
and  re-registering  of  ortho  photo  and  data  base  created  small  errors  in 
the  location,  orientation,  and  scale  of  any  given  registration  with  all 
others.  The  errorswere  not  apparent  to  the  eye,  but  when  radar  images 
were  simulated  from  this  data  base  and  the  cross  correlation  calculated 
between  these  simulated  and  real  images,  the  decrease  in  geometric  fidelity 
reduced  the  height  and  broadened  the  peak  of  the  cross  correlation.  This 
problem  also  had  the  potential  for  causing  a shift  in  the  location  of  the 
match  point,  the  peak  of  the  cross  correlation,  creating  an  apparent 
centering  error. 

This  problem  could  have  been  averted  easily  if  the  scale  of  the  ortho 
photo  had  been  reduced  photographically  to  that  desired  and  the  reduced 
ortho  photo  used  as  source  intelligence  for  the  data  base.  This  was  not 
done  because  the  table  was  available  and  past  experience  using  the  table 
for  other  kinds  of  projects  did  not  flag  a potential  problem.  In  all 
subsequent  work,  when  scale  reduction  has  been  desirable,  It  has  been 
realized  photographically. 
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It  was  impossible  to  fix  the  data  base  having  these  errors.  The 
only  correction  possible  would  have  been  scrapping  the  data  base  contain- 
ing the  problems  and  making  a new  one,  and  this  was  prohibitive  in  terms 
of  time  and  resources. 

6.2  Mis-Registration  of  Planimetry  and  Elevation  Data 

When  setting-up  the  coordinate  system  in  which  the  hand-drawn 
planimetry  for  the  bands  3 and  4 data  base  was  to  be  digitized,  one  point 
was  incorrectly  specified.  The  intent  was  to  digitize  the  planimetry 
data  in  the  coordinate  system  of  the  elevation  data.  This  was  desired 
to  minimize  problems  associated  with  merging  the  two  sets  of  data, 
planimetry  and  elevation.  Incorrect  specification  of  one  point  created 
both  scale  and  rotation  errors  of  the  digitized  planimetry  with  respect 
to  the  elevation  data. 

This  problem  was  not  detected  until  all  source  data  which  might 
have  contained  details  of  the  error  were  destroyed.  In  the  abscence  of 
facts  relating  to  the  exact  nature  of  the  error,  it  was  attempted  to 
affect  a scale  and  rotation  that  would  map  points  from  the  planimetry 
coordinate  system  into  the  desired  point  in  the  elevation  data  coordinate 
system.  A unique  mapping  was  impossible  to  obtain  because  of  the  problems 
with  the  geometric  fidelity  previously  discussed  (Section  6.1).  A best 
registration  between  planimetry  and  elevation  data  was  obtained  for 
selected  points  which  were  readily  identifiable  in  both  data  sets. 

Radar  simulations  were  created  from  the  data  base  containing  the 
merge  of  the  scaled  and  rotated  planimetry  with  the  elevation  data.  The 
cross  correlation  was  calculated  between  these  simulations  and  real  images. 
Misregistration  between  planimetry  and  elevation  data  caused  several 
problems  with  the  cross  correlation.  First,  the  location  of  the  match 
point  (the  peak  of  the  cross  correlation)  was  shifted  creating  an  apparent 
centering  error.  Second,  the  weight  was  reduced  and  the  width  broadened 
of  the  peak.  These  problems  mean,  of  course,  that  the  data  base  for  bands 
3 and  4 was  not  suitable  for  testing  the  PSM  as  a viable  radar  simulation 
technique  for  guidance  systems  using  the  Correlatron. 


7.0  RESULTS  AND  RECOMMENDATIONS 


7.1  Correlatron  Test  Results 

Radar  simulations  were  produced  for  four  (4)  different  scales  from 
the  two  (2)  data  bases  constructed  in  the  studies  reported  here.  Radar 
simulations  of  two  different  scales  were  made  from  each  of  the  data  bases. 
These  four  simulations  represented  the  altitudes  of  specific  points  on 
the  terminal  trajectory  of  a missile,  each  successively  lower.  The  highest 
altitude  simulation  is  called  band  4,  the  next  lower  is  called  band  3, 
next  is  band  2,  and  lowest  is  called  band  1.  Band  3 and  4 simulations 
were  produced  from  one  data  base,  and  band  1 and  2 simulations  from  the 
other. 

The  sequence  of  four  radar  simulations  were  tested  at  ETL  (U.S.  Army 
Engineer  Topographic  Laboratories)  by  a Correlatron.  The  simulations  were 
input  to  the  Correlatron  as  reference  guidance  scenes  for  correlation 
with  actual  radar  data  previously  collected  and  recorded  in  a flight  test 
program  conducted  over  the  Pickwick  test  site.  The  results  of  this  test 
are  mixed,  primarily  as  a result  of  the  problems  in  the  bands  3 and  4 
data  base  as  previously  discussed  (Section  6.0).  Complete  analyses  of 
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the  results  are  not  reported  here  as  that  is  beyond  the  scope  of  this 
report.  Qualitatively,  the  results  obtained  were: 

1)  Band  1 (lowest  altitude) 

Very  good  performance  within  accuracy  requirements; 

2)  Band  2 (higher  altitude) 

Excellent  performance  well  within  desired  accuracies; 

3)  Band  3 (next  higher  altitude) 

Failure  to  perform  within  accuracy  requirements; 

4)  Band  4 (highest  altitude) 

Acceptable  performance  within  accuracy  requirements. 

Band  1 and  2 simulations  were  developed  from  the  data  base  having  the 
highest  degree  of  detail.  Performance  of  these  simulations  was  spectacular 
in  view  of  the  fact  that  they  worked  this  well  on  the  first  trial. 


Simulations  produced  by  any  other  technique  have  required  numerous  trials, 
each  successive  trial  testing  a change,  to  obtain  acceptable  performance 
at  this  scale.  Simulations  produced  via  the  PSM  should  work  this  well 
as  the  PSM  rigorously  models  the  radar  guidance  system. 

Band  3 and  4 simulations  were  developed  from  the  data  base  having 
the  problems  discussed  in  Section  6.  Performance  of  these  simulations  was 
better  than  expected  in  view  of  the  nature  of  the  problems  in  the  data 
base.  The  problems  in  the  data  base  are  unfortunate  because  they  cloud 
what  otherwise  should  have  been  a clear  decision  in  favor  of  the  PSM. 
Acceptable  performance  has  been  readily  attained  by  most  other  simulation 
techniques,  in  part  because  of  the  large  reduction  in  detail  and  scale. 
Surely  if  the  band  3 and  4 data  base  had  not  contained  geometrical  errors, 
the  PSM  simulation  for  these  bands  would  have  performed  as  well  as  the 
band  1 and  2 simulations  did.  Of  course,  this  cannot  be  proved  with  the 
data  that  exist,  but  it  is  a reasonable  conjecture,  especially  in  view  of 
the  fact  that  even  with  the  errors, the  simulation  for  band  4 was  acceptable. 


7.2  Recommendations 


As  the  test  conducted  here  was  not  decisive,  it  is  recommended  that 
the  band  3 and  4 data  base  oe  rebuilt,  and  new  simulations  produced  and 
tested  on  the  Correlatron.  This  is  recommended  because  only  those  inti- 
mately familiar  with  the  technical  aspects  of  the  radar  simulation  problem 
will  recognize  that  the  band  3 failure  to  meet  expectations  was  for  the 
reasons  enumerated,  and  an  objective  test  needs  to  be  conducted  to  satisfy 
others.  This  is  important  because  the  PSM  offers  a technique  for  radar 
simulation  which  does  not  require  target  specific  special  processing  and 
fooling  around,  as  other  techniques  do.  All  the  PSM  requires  is  an  accurate 
data  base,  a requirement,  incidentally,  which  is  fundamental  to  any 
simulation  technique. 

It  is  further  recommended  that  the  PSM  be  used  to  produce  radar  simu- 
lations of  other  sites  for  testing  on  the  Correlatron.  Results  of  such 
tests  will  establish  the  validity  of  the  PSM  for  application  to  all  manner 
of  tarqet  scenes,  thereby  eliminating  the  site  specific  customizing  pre- 
**ntly  required  for  other  techniques. 
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Backscatter  Data 
for 

Bands  1 and  2 
Reference  Radar  Scenes: 
Pickwick  Site 


All  of  these  data  represent  the  following  radar  parameters: 
1)  Frequency  - X-band 


2)  Polarization  - Horizontal  transmit.  Horizontal  receive  (HH) 
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APPENDIX  B 


Backscatter  Data 
for 

Bands  3 and  4 
Reference  Radar  Scenes 
Pickwick  Site 


All  of  these  data  represent  the  following  radar  paramters: 

1)  Frequency  - X-band 

2)  Polarization  - Horizontal  Transmit,  Horizontal  receive  (HH) 
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POLAR  CONVERSION 

A.  POLAR  CREATE 

B.  POLAR  ARRAY 

C.  ARRAY  FIX 
REFERENCE  SCENE 

A.  POWER 

B.  ORE Y TONE 
RECTANGULAR  CONVERSION 

A.  RECTANGULAR  CREATE 

B.  RECTANGULAR  ARRAY 
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COMPUTER  PROGRAM 
I.  POLAR  CONVERSION 


This  computer  program  was  written  in  FORTRAN 
for  implementation  on  a Honeywell  66/60. 


It  consists  of 

three  subprograms 

A. 

Polar  Create 

B. 

Polar  Array 

C. 

Array  Fix 
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COMPUTER  PROGRAM 

II.  REFERENCE  SCENE 

This  computer  program  was  written  in  FORTRAN 
for  implementation  on  a Honeywell  66/60. 

It  consists  of  two  subprograms: 

A.  Power 

B.  Greytone 
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COMPUTER  PROGRAM 


III.  RECTANGULAR  CONVERSION 


This  computer  program  was  written  in  FORTRAN 
for  implementation  on  a Honeywell  66/60. 


It  consists  of  two  subprograms: 

A.  Rectangular  Create 

B.  Rectangular  Array 
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